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PART  I 


COMPOSITE  COMPUTER  PROGRAMS  FOR  EMPIRICAL 
EVALUATION  OF  INFRARED  SHIELD  CONFIGURATIONS 

SUMMARY 


Two  ray-tracing  programs  for  digital  computers  h»ve  been  written  to 
predict  the  enguler  distribution  of  radiant  intensity  from  an  infrared 
source  when  shielded  by  a  variety  of  cold  shields.  The  laboratory 
measurements  and  initial  prediction  efforts  were  reported  in  full  in 
the  Interim  report,  dated  October  1963,  of  this  contract. 

Studies  were  conducted  on  two  computers,  an  IBM  1620  and  an  IBM  7094. 
The  results  of  16  runs  on  the  1620  and  7  runs  on  the  7094  are  presented 
in  tabular  and  graphical  form  and  are  compared  with  the  laboratory 
measurements,  with  the  exception  of  the  last  four  experimental  runs, 
which  did  not  have  laboratory  counterparts.  These  studies  pointed  out 
the  power  of  this  technique  as  a  design  tool  and  Indicated  the  areas 
that  still  need  concentrated  work. 


CONCLUSIONS 


The  present  programs  did  not  solve  all  the  problems  associated  with  the 
ray-tracing  technique  of  determining  angular  radiation  distribution 
patternsi  but  improvements  were  made,  experience  was  gained,  and  the 
areas  needing  further  development  were  more  clearly  delineated.  The 
power  of  the  technique  as  a  design  tool  was  demonstrated  in  part  In  the 
last  four  runs. 

Further  study  is  needed  to  develop  a  suitable  correlation  factor  for 
relating  the  simple  techniques  of  ray  tracing  to  the  more  complex  physi¬ 
cal  processes  of  emission  and  reflection.  This  need  is  more  prominent 
at  large  angles. 

A  new  area  of  reeded  study  was  revealed,  namely,  the  detailed  dependence 
of  reflection  on  surface  characteristics  and  how  partially  diffuse  re¬ 
flections  can  be  handled  most  efficiently  in  a  ray-tracing  program. 

This  study  may  require  laboratory  investigations  to  establish  standards 
for  the  empirical  or  theoretical  treatment. 
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INTRODUCTION 


This  Is  the  final  report  on  a  study  of  the  measurement  and  prediction 
of  infrared  radiation  patterns  from  shielded  infrared  sources,  conducted 
by  Hayes  International  Corporation.  The  laboratory  measurements  pro¬ 
gram  and  the  Initial  efforts  at  pattern  prediction  by  a  digital  compu¬ 
ter  are  comprehensively  discussed  In  the  Interim  report*  and  will  not 
be  repeated  in  this  report.  Reference  should  be  made  to  the  interim 
report. 

That  portion  of  the  whole  program  which  Is  covered  by  this  report  arose 
from  certain  conclusions  drawn  In  the  Interim  report.  At  that  time, 
computer  programs  had  been  written  which  would  predict  radiation  patterns 
for  a  particular  set  of  shields  of  simple  geometry.  The  degree  of  agree¬ 
ment  with  laboratory  data  was  encouraging,  but  sufficient  discrepancies 
existed  for  some  of  the  configurations  to  justify  (1)  further  investi¬ 
gation  along  the  line  of  emoirical,  ray-tracing  computer  programs,  and 
(2)  development  of  an  analytical  formulation  using  integral  equations. 
Accordingly,  this  report  will  be  divided  into  two  parts,  the  first 
dealing  with  the  ray-tracing  program  and  the  second  with  the  integral 
equation  formulation. 


*  "A  Study  of  the  Radiative  Characteristics  of  Shielded  infrared 

Sources",  TREOOM  TR  63-45,  prepared  by  Hayes  International  Corp.  for 
the  U.  S.  Army  Transportation  Research  Comaand,  Fort  Eustis,  Virginia, 
October  1963, 
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EMPIRICAL.  RAY -TRACI KG  PKXMM6 


PROGRAM  LOGIC 


The  original  ray-tracing  programs,  discussed  in  tha  intar la  report, 
provided  for  ray  •  a  lection  in  tha  aperture  of  tha  model  and  geometrical 
tracing  back  to  tha  source  plate.  All  raye  ware  selected  ae  as  to  lie 
In  horizontal  planes.  Shan  the  new  computer  prograa*  were  written,  this 
original  approach  was  altered  so  that  the  progre*  would  more  closely 
parallel  the  phyaical  processes. 

The  origins  of  all  rays  were  selected  at  random  points  In  the  plane  of 
the  nodal  base,  corresponding  to  the  source  plate  of  the  laboratory 
apparatus.  To  each  ray  there  was  imparted  a  randoai  direction.  The  dis¬ 
tribution  of  randoai  directions  was  weighted  to  produce  e  cosine  distribu¬ 
tion  about  the  axis  of  ths  oodsl.  The  rays  were  net  restricted  to  horizontal 
planes.  As  before,  all  reflections  were  treated  as  specuhAr.  Generally,  a 
ray  was  traced  until  it  wea  reflected  back  to  the  base  or  out  through  the 
aperture.  Vhen  e  ray  passed  through  the  aperture  plane,  it  was  circularly 
projected  onto  the  horizontal  plane  by  converting  the  direction  cosine 
relative  to  the  axle  into  an  angle  relative  to  the  axis. 

Laboratory  measurements  had  indicated  that  the  radiating  source  plate 
approximated  a  Lambert  surface  within  an  average  6.2JJ  when  sampled  at  10° 
Intervals  from  0°  to  60°  and  at  65°.  Inspection  showed  that  the  surface 
emission  could  bo  better  approximated  by  a  quartlc  expression  than  by  the 
Lambertian  cosine  function.  An  expression  of  the  form  a/0*  +  b/r  +  c,  where 
jb  is  the  angle  between  the  rey  and  the  axis,  was  solved  for  suitable  values 
of  the  constants,  a,  b,  and  c.  A  plot  of  the  attenuating  function. 


f(j0)  =  .O7O104  -  .178^  +1,  (1) 


compared  with  the  laboratory  measurements  end  the  cosine  curve  is  shown 
in  Figure  1.  The  radiant  intensity  of  each  ray  was  attenuated  at  the  base 
plane  by  this  expression  according  to  its  sngle  of  emission. 

Eight  different  model  configurations  were  studied |  namely,  a  cylinder,  a 
cylinder  with  a  spherical  innerbody,  e  truncated  cone  with  the  aperture 
at  the  larger  end  (referred  to  as  the  large  cone),  the  large  cone  with  a 
spherical  innerbody,  a  truncated  cone  with  the  aperture  at  the  smaller  end 
(referred  to  as  the  smsll  cone),  the  small  cone  with  a  spherical  innerbody, 
a  model  composed  of  two  opposed  truncated  cones  having  their  large  diameters 
in  ceemon  (referred  to  ea  the  double  cone),  and  the  double  cone  with  a 
spherical  innerbody.  Two  programs  were  written*  one  fr r  the  cylindrical 
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Attenuation 
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and  simple  conical  models  with  and  without  innerbodies,  and  one  for  the 
double  cone  with  and  without  innerbody.  Both  programs  were  written  with 
provisions  for  changing  the  dimensions,  angles,  innerbof  size  and  posi¬ 
tion,  and  reflectivities  of  the  various  surfaces  so  that  a  variety  of 
configurations  could  be  examined.  Block  flow  diagrams  of  the  two  pro¬ 
grams  are  shown  in  Figures  2  and  3. 

The  confined  model  program,  written  for  the  cylinder  and  cones,  starts 
at  the  double-boxed  rectangle  in  the  upper  left  corner.  A  random  point 
is  located  in  the  base  plane  and  a  random  direction  is  imparted  to  an 
emitted  ray.  Tha  normal  radiant  intensity  is  attenuated  according  to 
the  angle  of  emission.  For  the  case  of  the  cylinder,  path  control  digits 
Lj,  and  1. 2  are  set  to  unity.  If  the  ray  intersects  the  innerbody  (ball), 
the  direction  cosines  of  the  reflected  ray  are  determined,  and  the  radiant 
intensity  is  attenuated  according  to  the  reflectivity  of  the  surface.  If 
the  direction  cosine  relative  to  the  axis  (m)  is  negative  and  the  model  Is 
the  cylinder  or  smell  cone,  a  new  direction  at  the  base  plane  is  assigned 
and  the  process  repeated.  If  m  is  positive  and  the  ray  does  not  pass 
through  the  exit  apeiture,  it  must  intersect  the  model.  In  the  case  of 
the  cylinder,  a  test  is  made  to  determine  how  fast  the  ray  is  progress¬ 
ing  down  the  model.  This  test  is  made  only  after  the  second  Intersec¬ 
tion  *ith  the  cylinder.  If  the  model  is  the  cylinder  or  1 f  the  m-direc- 
ticn  osine  is  positive,  a  test  for  intersection  with  the  innerbody 
made.  If  the  ray  doe*  not  intersect  the  innerbody  but  does  pass  through 
the  aperture,  the  emergent  ray  is  sorted  into  one  of  19  directions  accord¬ 
ing  to  its  angle  of  emergence.  When  a  predetermined  number  of  rays  have 
been  traced,  the  accumulated  data  is  output.  A  sample  output  date  sheet 
is  shown  in  Table  1.  PHI  is  the  angle  of  emergence.  Rays  emerging  at  a 
nominal  angle  +  2.5°  are  listed  under  the  nominal  angle.  For  example,  a 
ray  emerging  at  26.7°  would  be  Included  in  the  rays  at  25°,  while  a  ray 
at  27.7°  would  be  Included  in  the  rays  at  30°.  The  rays  at  0°  and  90 p 
have  only  a  2.5°  spread  Instead  of  the  usual  5°.  ARI  stands  for  average 
radiant  intensity  and  is  determined  by  dividing  TRI  (total  radiant  in¬ 
tensity)  by  NO.  RAYS.  The  line  at  the  bottom  of  the  print-out  gives  the 
specifications  for  the  particular  configuration  being  run.  The  abbrevia¬ 
tion  RIO  stands  for  the  initial  radiant  intensity  (at  the  base  plane ); 

RB,  for  radius  of  Innerbody  (ball)f  RC,  for  radius  of  cone  or  cylinder 
at  the  aperture  1  RP,  for  radius  of  cone  or  cylinder  at  the  source  plate; 
CB,  for  distance  along  r-axis  from  the  origin  to  the  center  of  the  Inner- 
body;  CLEN,  for  length  of  cone  or  cylinder;  REFB,  for  reflectivity  of 
cone  or  cylinder;  RAYS,  for  total  number  of  rays  traced  through  the  aper¬ 
ture;  and  ALFA,  for  the  naif-cone  angle  in  radians. 

The  double  cone  program  is  quite  similar  to  the  combined  model  program, 
except  tests  must  also  be  made  for  intersection  with  the  midplane  (the 
apex  of  the  double  cone),  located  at  a  distance  along  the  axis  of  the 
model  from  the  base  plane.  Cone  1  is  the  portion  of  the  model  between 
the  source  and  the  midplane,  and  Cone  2  is  the  portion  between  the  mid¬ 
plane  and  the  exit  aperture. 
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Figure  2.  Flow  Diagram:  Cylindrical  and 
Conical  Models. 
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Figure  3.  Flow  Diagram:  Double  Cone  Models 


TABLE  1 


SAMPLE  DATA  PAGE 


Radiation  Distribution  Program  -  Cy,  nder 
PHI  (DEG)  ARI  NO.  OF  RAYS  TRI 


0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

RIO  RB 

381819E-02  1.500 


. 3R1630E+02 
.345476E+02 
.  320656E+02 
.311343E+02 
.296948E+02 
.2  60420E+-02 
.235387E+02 
.2^3240E+02 
.202803E-02 
.203205E+02 
.  15344 9E+ 02 
.132982E+02 
.106233E+02 
.732909E+01 

.134207E+01 


RC  RP 

2.000  2.000 


220 

734 

1084 

1058 

1247 

873 

1084 

817 

611 

763 

324 

547 

273 

136 

136 

27 

0 

0 

0 

CB  CLEN 

6.3125  7.8125 


„ 839587E+04 
.253579E+05 
.347491E+05 
.329400E4-05 
. 370295E+05 
.227346E^05 
.255160E+05 
.182387E4-05 
.137297E+05 
.155045E+05 
.497174E+04 
.727412E+04 
. 29001 7E+04 
.  99675  7B-03 
.612055E+-03 
. 362360E+02 


REFB  REFC 

0.900  0.900 


RAYS  ALFA 

10000  0.0000 


RESULTS 


The  two  programs  were  originally  written  for  and  were  debugged  on  Hayes' 
IBM  1620  digital  computer.  Through  th  end  of  November  1963,  all  produc¬ 
tion  runs  were  also  run  on  the  1620.  By  that  time,  the  difficulty  of 
correlating  some  of  the  computer  predictions  with  the  laboratory  data 
suggested  the  need  for  a  faster  computer  in  order  to  trace  a  larger  num¬ 
ber  of  rays.  Consequently,  during  December  seven  runs  were  made  on  the 
IBM  7094  at  Eglin  Air  Force  Base,  Florida.  At  the  end  of  December,  five 
more  runs  were  made  on  the  1620. 

The  results  of  production  runs  made  through  November  are  shown  in  tabular 
form  in  Table  2.  The  corresponding  graphs  are  shown  in  Figures  4  through 
13,  in  which  the  laboratory  curves  were  taken  from  the  interim  report. 

The  number  of  rays  per  angle  interval  is  shown  as  an  indication  of  the 
reliability  of  the  data.  Data  averaged  over  20  to  30  rays  would  be  ex¬ 
pected  to  exhibit  less  scatter  than  data  averaged  over  3  or  4  rays. 

There  are  a  few  exceptions,  but  this  is  the  general  rule.  The  notation 
E+-XX  in  the  data  fields  is  a  computer  convention  for  powers  of  tens;  thus 
.152662E-01  stands  for  .152662  x  10"1. 

FigMre  4  shows  excellent  agreement  between  the  computer  prediction  and 
the  laboratory  measurement,  even  for  the  relatively  small  number  of  rays 
traced. 

The  agreement  in  Figure  5  is  better  than  it  appears,  because  the  discre¬ 
pancy  along  the  0°-line  is  somewhat  arbitrary.  The  quantity  RIO  In  the 
computer  program  represents  the  initial  radiant  intensity  along  the  axis 
of  the  cylinder  and  is  input  data  to  the  programs.  The  actual  value  used 
in  most  cases  was  obtained  by  averaging  a  large  number  of  laboratory  meas¬ 
urements.  The  value  for  a  particular  laboratory  run  usually  varied  some¬ 
what  from  this  average. 

When  normalized  to  the  laboratory  measurement,  the  prediction  of  Figure  6 
agrees  around  to  40°,  diverging  thereafter  to  a  difference  at  70°  by  a 
factor  of  five.  Part  of  this  discrepancy  may  be  ittributable  to  inac¬ 
curacy  in  details  of  the  quart! c  attenuation  function  at  large  angles, 
since  these  details  had  to  be  determined  by  extrapolation  for  angles 
greater  than  60°.  Furthermore,  at  all  angles  over  50°,  no  more  than  two 
rays  were  traced.  The  results  of  so  few  tracings  can  hardly  be  expected 
to  give  more  than  a  rough  estimate.  This  latter  problem  was  largely 
overcome  by  going  to  the  7094  computer.  Finally,  the  presence  of  an 
innerbody  tends  to  attenuate  the  radiation  more  than  simple  ray-tracing 
predicts . 

The  configuration  represented  by  Figure  7  has  been  the  subiect  of  con¬ 
centrated  efforts,  partly  because  it  is  the  fastest  program  and  partly 
because  of  the  sharp  cut-off  characteristic  near  40°.  One  of  the  big 
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TABLE  2 


RADIATION  DISTRIBUTION  DATA 
(IBM  1620  RUNS) 


Configuration 

and 

Specification 

0 

(degrees ) 

Averaqe 

Radiant  Intensity 
(watts/»ter.) 

No.  of  Ray* 

Gold  Cylinder 

0 

.381657E+02 

4 

RIO  .3818I9E-02 

5 

.380185E+02 

17 

RB  0 

10 

. 373501 E+02 

13 

RC  2.0 

15 

.361472E+02 

12 

RP  2.0 

20 

. 344524E+02 

0 

CB  0 

25 

.331676E+02 

4 

CLEN  7.8125 

30 

.314157E+02 

15 

RHFB  0 

35 

.2889388+02 

10 

REFC  .969 

40 

.267623E+02 

20 

RAYS  200 

45 

.238183E+02 

9 

ALFA  0 

50 

.217S16E+02 

12 

55 

.18359* +02 

5 

60 

.1584O7E+02 

6 

65 

. 12664 3E +02 

12 

70 

.951432E+02 

10 

75 

•  659243E+02 

16 

80 

.3537361+02 

10 

85 

.  145087E+01 

11 

90 

.239796E+00 

6 

Black  Cylinder 

0 

.381656E+02 

15 

RIO  .  38181 9E- 02 

5 

.276907E+02 

33 

RB  0 

10 

.279422E+02 

31 

RC  2.0 

15 

.176851E+02 

23 

RP  2.0 

20 

.492566E+01 

30 

CB  0 

25 

•  1710ME+00 

19 

CLEN  7.8125 

30 

.  126899E+00 

37 

REFB  0 

35 

. 987140E-  01 

23 

REFC  .005 

40 

.571691E-  01 

40 

RAYS  300 

45 

.641878E-03 

29 

ALFA  0 

50 

.576748E-03 

20 

II 


TABLE  2  (OONT’D.) 


RADI  AT  iCV  DISTRIBUTION  DATA 
(IBM  1620  RUNS) 


Configuration  .  Average 

and  “  Radiant  Intensity  No.  of  Raya 

Spaclficationa _ (dagraaa)  (itf/iter.) _ 


Gold  Cylinder  and 

0 

.3603048+02 

4 

Gold  Innerbody 

5 

.3534886+02 

3 

RIO 

.3818196-02 

10 

.338507E+02 

12 

RB 

1.5 

15 

.319862E+02 

17 

RC 

2.0 

20 

.  31581 4E+0T 

19 

RP 

2.0 

25 

.27O390E+02 

6 

CB 

6.3125 

30 

.269472E+02 

9 

CLEN 

7.8125 

35 

.2487246+02 

8 

REFB 

.950 

40 

. 188747E  +02 

5 

REFC 

.950 

45 

.226567E+02 

5 

RAYS 

100 

50 

.  1 63547E+02 

2 

ALFA 

0 

» 

.253561E+-02 

2 

60 

.1155866+  02 

1 

65 

.9264286+- 01 

2 

70 

.713797E+01 

1 

75 

.366703E+01 

2 

80 

.2232056+01 

2 

Gold 

Large  Cone 

0 

.588903E+02 

14 

RIO 

.381819E-02 

5 

.5629166+02 

45 

RB 

0 

10 

.5461496+02 

51 

RC 

2.5 

15 

.4659356+02 

72 

RP 

1.375 

20 

.3871828+02 

69 

C8 

0 

25 

.2503528+02 

90 

CLEN 

8.0 

30 

.1809868+02 

101 

REFB 

0 

3*. 

.148748E+02 

42 

REFC 

.95 

40 

.4381 72E+  01 

16 

RAYS 

500 

ALFA 

.1396 
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TABLE  2  (OOfT'D.) 


RADIATION  DISTRIBUTION  DATA 
(IBM  1620  RUNS) 


Configuration 

and 

Specifications 

0 

(degrees ) 

Average 

Radiant  Intensity 
(watts/ster.) 

No.  of  Rays 

Black  Large  Cone 

0 

.434506E+02 

22 

RIO 

.381819E-02 

5 

.269288E+02 

71 

RB 

0 

10 

. 303066E+02 

66 

RC 

2.5 

15 

.136240E+02 

85 

RP 

1.375 

20 

.342977E+01 

101 

CB 

0 

25 

.124669E+00 

85 

CLEN 

8.0 

30 

.637961E-01 

57 

REFB 

0 

35 

.125760E+00 

13 

REFC 

.006 

RAYS 

500 

ALFA 

.1396 

Gold 

Large  Cone  and 

0 

.40703CE+02 

3 

Gold  Innerbody 

5 

.495897E+Q2 

26 

RIO 

•35CE-02 

10 

.44923CE+02 

32 

RB 

1.75 

15 

.427347E+-02 

29 

RC 

2.5 

20 

. 308968E+02 

29 

RP 

1.375 

25 

.275685E+02 

32 

CB 

6.25 

30 

.21722 OE +02 

28 

CLEN 

8.0 

35 

.205582E+02 

8 

REFB 

.95 

40 

.258950E+02 

9 

REFC 

.95 

45 

.256665E+02 

4 

RAYS 

200 

ALFA 

.1396 

Gold 

Snail  Cone 

0 

.  180428E+C2 

22 

RIO 

.381819E-02 

5 

.179811E+02 

31 

RB 

0 

10 

.  177624E+02 

31 

RC 

1.375 

15 

.1718381+02 

32 

RP 

2.5 

20 

.169863E+02 

31 

CB 

0 

25 

.167545E+02 

32 

CLEN 

8.0 

30 

•165136E+02 

41 

REFB 

0 

35 

.160063E+02 

33 

REFC 

.95 

40 

.152817E+02 

24 

RAYS 

400 

45 

.152970B+02 

10 

ALFA 

-.1396 

50 

.148706E+02 

17 

13 


TABLE  2  (OOWT'D.) 
RADIATION  DISTRIBUTION  DATA 
(IBM  1620  RUNS) 


Configuration 

and 

Specifications 

0 

(degrees ) 

Average 

Radiant  Intensity 
(watts/ster. ) 

No.  of 

Gold  Snail  Cone 

55 

. 147235E+02 

14 

(Cont 'd. ) 

60 

.  134366E+02 

34 

65 

.132366E+02 

23 

70 

.  128122E+02 

19 

75 

- 

0 

80 

. 117826E+02 

6 

Gold  Small  Cone  and 

0 

.  180433E+02 

11 

Gold  Innerbody 

5 

.177131E+02 

20 

RIO  .381819E-02 

10 

.173296E+02 

16 

RB  1.0 

15 

. 161585E+02 

33 

RC  1.375 

20 

.162184E+02 

42 

RP  2.500 

25 

.159515E+02 

30 

CB  7.0 

30 

.156533E+02 

47 

CLEN  8.0 

35 

.147551E+02 

28 

REFB  .95 

40 

. 149061E+O2 

19 

REFC  .95 

45 

.148379E+02 

28 

RAYS  400 

50 

.139076E+02 

33 

ALFA  -.1396 

55 

. 137101E+02 

22 

60 

.133707E-+02 

30 

65 

.  131336E-+02 

11 

70 

•  121432E-HX2 

6 

75 

.113070E-K32 

12 

80 

.109579E+02 

9 

85 

.845597E-K31 

3 

Gold  Double  Cone 

0 

.313432E+02 

10 

RIO  .381819E-02 

5 

.304695E+02 

46 

RC  1.8125 

10 

.290311E+02 

23 

RM  2.3437 

15 

. 28001 4E+02 

17 

RP  1.5625 

20 

.225717E-KJ2 

30 

RB  0 

25 

.  198778E-KJ2 

23 

CB  0 

30 

.177679E+02 

36 

CLEN  5.5 

35 

.248925E-KJ2 

26 

ZM  2.9375 

40 

.217852E+02 

29 

REFC  1  .95 

43 

.247496E+02 

29 

REFC  2  .95 

50 

.202480E+02 

40 

REFB  0 

55 

.206132E+02 

41 

14 


TABLE  2  (CONT  'D. ) 
RADIATION  DISTRIBUTION  DATA 
(IBM  1620  RUNS) 


Configuration 

and 

Spec! fications 

0 

(degrees ) 

Average 

Radiant  Intensity 
( watts/s ter. ) 

No.  of  Rays 

Geld  Double  Cone 

60 

.168658E+02 

16 

(Cont 'd. ) 

65 

.107599B-02 

11 

RAYS  400 

70 

. 122692E+02 

4 

ALFA  1  .2598 

75 

.  13954  0E+02 

3 

ALFA  2  -.2042 

80 

.  1 10883E+02 

8 

85 

.174685E+02 

3 

Gold  Double  Cone  and 

0 

_ 

C 

Gold  Innerbody 

5 

.187669E+02 

13 

RIO  . 381819E-02 

10 

.207801E+02 

22 

RC  1.8125 

15 

. 190052E+02 

14 

RM  2.3437 

20 

. 198098 B-  02 

27 

RP  1.5625 

25 

.216940E+02 

38 

RB  1.75 

30 

•207317E+02 

53 

CB  2.9375 

35 

.190096E+02 

49 

CLEN  5.5 

40 

.219302E+02 

41 

ZM  2.9375 

45 

.201971E+02 

35 

REFC  1  .95 

50 

.  18282 0E+ 02 

41 

REFC  2  .95 

55 

.150247E+02 

35 

REFfc  .95 

60 

.144817E+02 

13 

RAYS  400 

65 

.117407E+02 

12 

ALFA  1  .2598 

70 

.300514E+C1 

4 

ALFA  2  -.2042 

75 

.148839E+02 

1 

80 

- 

o 

85 

.  123987E+02 

2 

15 
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Figure  ^ .  Black  Large  Cone  without  Figure  9.  Gold  Large  Cone  with  Gold 

Innerbody.  Innerbody. 


lOftO 
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Figure  10.  Gol  .1  Small  Con**  *if  i.tu'  Fiquro  11.  Gol  i  Small  Ccne  wit!  Gol 

I r.n.»»rbo  iy  .  I-v.erbo.iy. 
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Figure  12.  Gold  Double  Cone  without  Figure  13.  Gold  Double  Cone  with  Gold 

Innerbody.  Innerbody. 


problems  has  been  the  slight  bulge  between  30°  and  35°,  which  becomes 
quite  pronounced  when  the  innerbody  is  added.  A  slight  change  in  the 
program,  to  be  described  presently,  reduced  this  problem  somewhat. 

Figure  8  shows  the  effect  of  a  low-reflectivity  surface  material.  Except 
for  the  on-axis  discrepancy,  the  two  curves  probably  could  be  brought  into 
good  agreement  by  varying  the  value  of  the  surface  reflectivity. 

Figure  9  represents  the  large  cone  with  innerbody.  As  mentioned  above, 
the  hump  at  40°  to  45°  is  the  biggest  problem.  No  completely  satisfactory 
solution  to  this  problem  has  been  devised. 

Figures  7  and  10  played  a  prominent  part  in  the  formulation  of  an  early 
empirical  correction  factor,  designed  to  correct  the  diverging  discrepancy 
at  large  angles.  The  correction  factor  was  of  the  form 

f(0>  -j-f-)  =  [cos  A 

in  which  the  quantities  0,  RC,  and  RP  have  the  same  significance  as  in 
the  tabulated  data.  This  multiplicative  factor  had  the  effect  of  boost¬ 
ing,  with  increasing  0,  the  radiant  intensity  for  the  large  cone  and  of 
attenuating  that  of  the  small  cone,  while  leaving  the  cylinder  unchanged. 

The  factor,  however,  was  not  applicable  when  innerbodies  were  present  or 
for  double  cone  configurations.  Extensive  efforts  were  made  to  develop 
an  empirical  expression  which  would  correlate  predictions  with  laboratory 
measurement  for  all  configurations  simply  by  substituting  the  appropriate 
values  of  certain  geometrical  parameters.  These  attempts  were  unsuccessful. 

The  characteristic  feature  of  the  laboratory  curve  for  the  double  cone 
(Figure  12)  is  the  depression  at  40°.  This  configuration  has  been  run 
several  times,  and  the  depression  always  occurs  somewhere  between  25°  and 
30°.  It  is  felt  that  the  points  beyond  65°  are  not  reliable,  because  the 
number  of  rays  traced  dropped  sharply  in  this  region. 

When  normalized  to  the  laboratory  on-axis  value,  the  prediction  shown  in 
Figure  13  agrees  well  with  the  experimental  curve  around  to  about  50°. 

Again,  the  points  beyond  65°  are  susceptible  to  suspicion  because  of  the 
paucity  of  rays  traced. 

The  problem  of  correlating  prediction  curves  with  laboratory  curves  led 
to  repeated  efforts  to  develop  a  suitable  correlation  factor.  It  did 
not  seem  unreasonable  that  such  a  factor  might  be  needed,  for  a  simple 
macroscopic  averaging  technique  was  being  used  to  simulate  a  complex 
quantum  process.  It  was  thought  that  perhaps  the  distribution  of  rays 
might  provide  a  statistical  correction  factor,  but  to  explore  this 
possibility,  it  was  necessary  to  trace  many  times  the  number  of  rays 
that  could  be  traced  on  the  1620  computer.  It  has  been  noted,  also, 
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that  the  small  number  of  rays  traced  might  be  responsible  for  the  large 
degree  of  scatter,  particularly  at  large  angles.  In  order  to  trace  the 
large  number  of  rays  needed  In  economically  feasible  times,  it  was  neces¬ 
sary  to  go  to  a  faster  computer.  For  this  task,  the  IBM  7094  at  Eglin 
Air  Force  Base,  Florida,  was  chosen.  Seven  configurations  were  run,  each 
for  10,000  rays.  It  turned  out  that  one  run,  which  would  have  required 
250  hours  on  the  1620,  was  completed  on  the  7094  in  less  than  7  minutes. 
Another  run,  requiring  60  hours  on  the  1620,  ran  to  completion  in  2.1 
minutes  on  the  7094.  The  tabulated  data  from  these  runs  is  shown  In 
Table  3.  The  corresponding  graphs  are  shown  in  Figures  14  through  20. 

Figure  14  bears  a  close  resemblance  to  Figure  4,  except  the  scatter  is 
reduced  and  the  tail  of  the  prediction  curve  is  about  65#  closer  to  the 
laboratory  curve.  This  points  up  the  magnitude  of  change  that  is  possible 
when  a  relatively  large  number  of  rays  are  traced. 

In  Figure  15  the  agreement  is  improved  at  30°,  but  worsened  at  40°. 

Figure  16  shows  an  improvement  over  the  corresponding  earlier  graph, 

Figure  8.  When  normalized  to  the  on-axis  laboratory  value,  the  two 
curves  of  Figure  14  nearly  coincide. 

Figures  17  and  18  are  essentially  the  same  as  the  earlier  predictions. 

The  main  difference  in  Figure  19  Is  in  the  region  of  65°  to  70°.  The 
points  beyond  70°  were  inexplicably  high.  Finally,  Figure  20  was  essen¬ 
tially  unchanged. 

After  the  completion  of  these  runs,  a  change  was  made  in  the  base  plane 
attenuation  function  (Eq.  1).  Figure  21  shows  the  quartic  attenuation 
function  of  Eq.  1  and  the  estimated  laboratory  curve.  Since  a  number  of 
the  predictions  were  unexpectedly  high  at  large  angles  (greater  than  70°), 
it  w'.s  felt  that  altering  the  attenuation  function  to  better  fit  the 
laboratory  curve  might  improve  the  prediction  curves.  Two  terms  were 
added,  giving  the  attenuation  function  the  form 

f(0)  =  -.01463  06  +  .0701  04  -.5786  ft  -.22  sin10  0+1.  (2) 


The  revised  function  is  also  shown  in  Figure  21.  Two  short  runs  were 
made  on  the  large  cone  configurations  (Figures  22  and  23).  The  con¬ 
figuration  without  innerbody  showed  a  slight  improvement  at  40°,  but 
the  configuration  with  innerbody  (158  rays  traced)  showed  so  much  scatter 
that  a  conclusive  curve  could  not  be  drawn.  Notably,  however,  the  range 
was  extended  to  50°. 

At  the  conclusion  of  these  studies,  four  additional  runs  were  made  on 
the  double  cone  program,  simulating  configurations  for  which  there  were 
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TABLE  3 


RADIATION  DISTRIBUTION  DATA 

(IBM  7094  RUNS) 

Configuration 

H 

Average 

and 

y> 

Radiant  Intensity 

No.  of  Rays 

Specifications 

(degrees ) 

(watts/ster. ) 

Gold  Cylinder  and 

0 

. 381630E+02 

220 

Gold  Innerbody 

5 

.345476E+02 

734 

RIO 

.381819E-02 

10 

.320656E+02 

1084 

RB 

1.5 

15 

.311343B-02 

1058 

RC 

2.0 

20 

.296948&4-02 

1247 

RP 

2.0 

25 

.260420E+-02 

873 

CB 

6.3125 

30 

.235387E+02 

1084 

CLEN 

7.8125 

35 

.223240&-02 

817 

REFB 

.90 

40 

•  202803E+Q2 

677 

REFC 

.90 

45 

.203205E+02 

763 

RAYS 

10000 

50 

.153449&-02 

324 

ALFA 

0 

55 

.132982B-02 

547 

60 

.106233E+02 

273 

65 

. 732909 E+-01 

136 

70 

.450041&+-01 

136 

75 

.134207E+01 

27 

Gold 

Large  Cone 

0 

.580217E+02 

352 

RIO 

.381819E-02 

5 

.563407E+02 

1564 

RB 

0 

10 

.536843E+02 

1732 

RC 

2.5 

15 

.461120E+02 

1647 

RP 

1.375 

20 

.403341E+02 

1500 

CB 

0 

25 

.310808E+02 

1489 

CLEN 

8.0 

30 

.225956E+02 

1026 

REFB 

0 

35 

.144971E-C02 

606 

REFC 

.90 

43 

.965673E+01 

84 

RAYS 

10000 

ALFA 

.1396 
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TABLE  3  (CONT ’D. ) 
RADIATION  DISTRIBUTION  DATA 
(IBM  7094  R'JNS ) 


Configuration 

and 

Specifications 

(degrees ) 

Average 

Radiant  Intensity 
(watts/ster. ) 

No.  of  Rays 

Black  Large  Cone 

0 

.436898E+02 

402 

RIO 

•381819E-02 

5 

.324186E+02 

1469 

RB 

0 

10 

.  30368 7E+ 02 

1921 

RC 

2.5 

15 

.  189675E+02 

1939 

RP 

1.375 

20 

.739769E+01 

1973 

CB 

0 

25 

.227655E+01 

1401 

CLEN 

8.0 

30 

.772978E+01 

713 

REFB 

0 

35 

.731115E+01 

182 

REFC 

.50 

RAYS 

10000 

ALFA 

.1396 

Gold  Large  Cond  and 

0 

. 40901 7E +02 

95 

Gold  Innerbody 

5 

. 49071 4E +02 

1235 

RIO 

.381819E-02 

10 

.451746E+02 

1666 

RB 

1.75 

15 

.382927E+02 

1632 

RC 

2.50 

20 

.313438E+02 

1527 

RP 

1.375 

25 

.252095E+02 

1426 

CB 

6.25 

30 

.185244E+02 

1290 

CLEN 

8.0 

35 

.157011E+02 

447 

REFB 

.90 

40 

. 1 67201 E +02 

517 

REFC 

.90 

45 

. 17899 OE +02 

165 

RAYS 

10000 

ALFA 

.1396 

Gold  Snail  Cone  and 

0 

.180414E+02 

286 

Gold  Innerbody 

5 

.177788E+02 

323 

RIO 

.381819E-02 

10 

.171025E+02 

467 

RB 

1.0 

15 

. 1 65894E+  02 

935 

RC 

1.375 

20 

. 163068 E+ 02 

1222 

RP 

2.5 

25 

.161458E+02 

648 

CB 

7.0 

30 

.155478E+02 

1050 

CLEN 

8.0 

35 

.152107E+02 

969 

REFB 

.95 

40 

.  150018E+02 

496 

REFC 

.95 

45 

.147765E+02 

791 

RAYS 

10000 

50 

.13976CE+02 

611 

ALFA 

-.1396 

55 

.  138936E+02 

728 

24 


TABLE  3  (CONT'D.) 
RADIATION  DISTRIBUTION  DATA 


(IBM  7094  RUNS) 

Configuration 

Average 

and 

(degrees ) 

Radiant  Intensity 

No.  of  Rays 

Specifications 

(watts/ster. ) 

Gold  Saall  Cone  and 

60 

.133401E+02 

504 

Gold  Innerbody  (Cont’d 

.)  65 

.128702E-KJ2 

395 

70 

.116750E+02 

106 

75 

. 112730E+02 

180 

80 

.108286E+02 

215 

85 

.908789E+01 

72 

Gold  Double  Cone 

0 

.3134868+02 

357 

RIO  .381819E-02 

5 

. 303865E+02 

920 

RC  1.8125 

10 

.3044206+02 

732 

RM  2.3437 

13 

.262457E+02 

578 

RP  1.5625 

20 

.239363E-KJ2 

682 

RB  0 

25 

. 193568E+02 

850 

CB  0 

30 

• 189916E+02 

919 

CLEN  5.5 

35 

.225223E+02 

749 

ZM  2.9397 

40 

.223933E+02 

787 

REFC  1  .95 

45 

.242015E+02 

629 

REFC  2  .95 

50 

.225313E+02 

834 

REFB  0 

55 

.209491E+02 

650 

RAYS  10000 

60 

. 175506E+02 

460 

ALFA  1  .2598 

65 

.126358E+02 

324 

ALFA  2  -.2042 

70 

.904595E+01 

119 

75 

.405550E+01 

120 

80 

.117768E+02 

137 

85 

.  196587E+02 

119 

90 

.139184E+02 

34 

Gold  Double  Cone  and 

0 

.246849E+02 

76 

Gold  Ball 

5 

. 21402 7E +02 

257 

RIO  .381819E-02 

10 

.  2  05 852 E  +02 

439 

RC  1.8125 

15 

.  19704 3E -KJ2 

567 

RM  2.3437 

20 

.  2 06 064 E +02 

849 

RP  1.5625 

25 

.223031E+02 

772 

RB  1.75 

30 

.204937E+02 

1005 

CB  2.9397 

3r 

.202852E  +02 

902 

CLEN  5 « 5 

40 

•215292E+02 

1333 

ZM  2.9397 

45 

.197642E+02 

1173 

25 


TABLE  3  (OONT 'D. ) 
RADIATION  DISTRIBUTION  DATA 
(IBM  7094  RUNS) 


Configuration 

and 

Specifications 

0 

(degrees ) 

Average 

Radiant  Intensity 
(watts/ster. ) 

No.  of  Rays 

Gold  Double 

Cone  and 

50 

.193257E+02 

1080 

Gold  Ball 

(Cent 'd. 

55 

.1453375+02 

645 

R5FC  1 

.9? 

60 

.1431168+02 

4  38 

RHFC  2 

.9b 

(■ 5 

.1187808+02 

310 

REFB 

.95 

70 

.5031 35 E+ 01 

77 

RAYS 

10000 

"6 

.1488405+02 

26 

ALFA  1 

.2598 

S  0 

- 

0 

ALFA  2 

-.2042 

85 

. 123988E+02 

26 

90 

.138592E+02 

25 

26 
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Gold  Small  Cc'to  with  Gold 
Innerbody 
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no  laboratory  counterparts.  These  configurations  ware  the  possible  com¬ 
binations  of  atodels  and  innerbodies  having  reflectivities  of  0.1  and  0.4. 

Figure  24  represents  the  configuration  with  both  aodel  and  innerbody 
having  reflectivities  of  0.1.  The  solid  curve  connects  all  points,  while 
the  broken  curve  suggests  an  alternative  plot  which  seeks  to  average  out 
rapid  fluctuations. 

Curve  B  in  Figure  25a  represents  the  configuration  in  which  the  conical 
sections  have  a  reflectivity  of  0.1  and  the  innerbody  has  a  reflectivity 
of  0.4.  The  broken  line  bears  a  strong  resemblance  to  the  broken  line 
of  Figure  24,  suggesting  that  the  reflectivity  of  the  innerbody  does  not 
seriously  affect  the  over-all  radiation  pattern.  For  curve  C  the  reflec¬ 
tivities  of  model  and  Innerbody  mere  interchanged,  and  the  curve  shows  a 
corresponding  upward  displacement.  Finally,  for  Figure  Jsb  all  reflec¬ 
tivities  are  0.4.  Around  to  50°  this  curve  closely  resembles  curve  C. 

The  points  around  70°  are  suspicious  because  only  single  rays  traced. 

For  all  these  four  configurations,  the  total  number  of  rays  traced  was 
small,  and  the  curves  reflect  the  expectedly  wide  scatter.  The  tabulated 
data  of  these  runs  is  presented  for  reference  in  Table  4. 
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TABLE  4 


RADIATION  DISTRIBUTION  DATA 
ADDITIONAL  DOUBLE  OONE  CONFIGURATIONS 


Configuration 

and 

Specifications 

$ 

(degrees ) 

Average 

Radiant  Intensity 
(watts/ster. ) 

No.  of  Rays 

Double  Cone  and 

0 

.299860E-03 

1 

Innerbody,  "A” 

5 

. 126905 E+01 

4 

RIO  .381819E-02 

10 

.987272E+00 

10 

PC  1.8125 

15 

.149825E*01 

18 

RM  2.3437 

20 

.103791E+01 

18 

RP  1.5625 

25 

. 453564 E+00 

34 

RB  1.75 

30 

.267939E+00 

27 

CB  2.9397 

35 

.349939B-00 

25 

CLEN  5.5 

40 

.162114E+01 

28 

ZM  2.9347 

45 

.  118030E+01 

26 

REFC  1  .10 

50 

.840194E+01 

23 

REFC  2  .10 

55 

.136106E-Q2 

8 

REFB  .10 

60 

.157073E-01 

4 

RAYS  232 

65 

. 106748E-01 

2 

ALFA  1  .2598 

70 

.207747E-01 

] 

ALFA  2  -.2042 

75 

.110112E-01 

3 

Double  Cone  and 

0 

.258076E+00 

2 

Innerbody,  "B" 

5 

.149914E+01 

5 

REFC  1  .10 

10 

.133406E+01 

8 

REFC  2  .10 

15 

.151143E+01 

17 

REFB  .40 

20 

.912819E+00 

22 

25 

.774297E+00 

13 

30 

.502498E+00 

19 

f 

35 

.504815E+00 

16 

40 

.  150352E+-01 

27 

45 

. 930338 E+ 00 

27 

50 

.863313E+00 

23 

55 

•381679E-01 

12 

60 

.197174E-01 

8 

65 

.330530E-01 

3 

TABLE  4  (OONT 'D. ) 


RADIATION  DISTRIBUTION  DATA 
ADDITIONAL  DOUBLE  CONE  CONFIOJRATIONS 


Configuration 

and 

Speci  f  ications 

0 

(decrees ) 

Average 

Radiant  Intensity 
(watts/ster. ) 

No.  of  Rays 

Double  Cone  and 

0 

.501221E-02 

1 

Innerbody,  "C" 

5 

. 601 125E+01 

5 

REFC  1  .40 

10 

.278228E+01 

9 

REFC  2  .40 

15 

•542002E+01 

17 

REFB  .10 

20 

.3193705*01 

23 

25 

.296846E+01 

20 

30 

.325373E+01 

22 

35 

.274256E+01 

27 

40 

.517727H+01 

30 

45 

. 337503E+01 

23 

50 

. 394718E+01 

28 

55 

.414241E+00 

15 

60 

.969764E-01 

8 

65 

.163193E+00 

w/ 

Double  Cone  and 

0 

. 320781E+00 

1 

Innerbody,  ”D" 

5 

. 695572E+01 

6 

REFC  1  .40 

10 

.283313E+01 

11 

REFC  2  .40 

15 

.509215E+01 

18 

REFB  .40 

20 

.414264 E+01 

22 

25 

.435535E+01 

18 

30 

.372848E+01 

25 

35 

. 304150E+01 

24 

40 

.492055E+01 

24 

45 

.400914E+01 

21 

50 

.407132E+01 

26 

55 

.516437E+00 

14 

60 

.353947E+00 

11 

65 

.547787E+00 

7 

70 

- 

0 

75 

.  1 1 1 103E+01 

1 

80 

- 

0 

85 

.690870E-01 

1 

90 

. 106347E-01 

1 

34 


soae 


3*3 


Figure  24.  Double  Cone  witn  Innerbody,  Figure  2c>a  .  Double  Cone  with  Innerbody 

Cur/e  A.  Curves  p  and  C. 


Radiant  Intensity  (watts/ater. ) 


10 


20 


30 


Figure  25b.  Double  Cone  with  Innerbody, 
Curve  D. 
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EVALUATION 


This  study  was  undertaken  in  an  attempt  to  improve  the  prediction  programs 
over  their  earlier  counterparts.  The  programs  were  made  as  flexible  as 
possible  without  sacrificing  unduly  their  operating  efficiency.  The  ray 
tracing  was  made  to  proceed  from  the  base  plane  outward,  since  this  paral¬ 
leled  the  physical  process.  The  random  numbers  used  to  locate  origins  of 
rays  and  to  determine  directions  were  taken  from  a  table  of  uniformly 
random  numbers,  but  only  2400  different  numbers  w«re  used.  The  list  was 
used  repeatedly  until  the  required  number  of  rays  were  traced.  A  more 
advantageous  arrangement  would  have  employed  a  random  number  generator. 

This  problem  did  not  become  significant  until  the  runs  were  made  on  the 
7094,  at  which  time  a  suitable  generator  was  not  readily  available. 

Examination  of  the  plots  revealed  that  the  programs  yielded  accurate  pre¬ 
dictions  for  some  configurations,  but  failed  to  give  completely  satisfac¬ 
tory  results  for  others.  Most  of  the  runs,  however,  fell  somewhere  in 

the  broad  region  between  these  two  extremes.  In  every  case,  there  was  at 
least  a  measure  of  qualitative  agreement. 

From  the  experience  gained  during  the  course  of  this  study  a  number  of 

areas  of  improvement  have  become  apparent.  The  most  urgent  need  is  for 

a  computer  faster  than  the  1620.  In  many  cases  the  discrepancies  might 
be  resolved  by  trying  various  combinations  of  parameters,  such  as  the 
initial  radiant  intensity  and  the  reflectivities.  Another  area  deserv¬ 
ing  study  is  the  treatment  of  reflections  from  surfaces  which  are  neither 
purely  diffuse  nor  purely  specular.  This  may  be  a  difficult  problem  and 
might  require  some  experimental  laboratory  work  to  provide  standards 
against  which  the  empirical  or  theoretical  models  could  be  judged.  There 
is  still  considerable  work  to  be  done  on  the  development  of  a  suitable 
correlation  factor  needed  to  compensate  for  fundamental  differences  be¬ 
tween  the  physical  processes  and  the  simplified  techniques  of  the  computer 
programs . 

In  spite  of  the  discrepancies  still  present  in  this  development,  there  is 
strong  evidence  that  this  approach  to  the  determination  of  radiant  in¬ 
tensity  patterns  can  be  developed  into  a  valuable  design  tool.  The  work 
accomplished  so  far  has  revealed  the  areas  requiring  more  study,  but  more 
important,  it  has  given  a  glimpse  of  the  power  of  the  tool.  Design  changes 
can  be  made  and  investigated  in  minutes  or  hours  instead  of  the  days  or 
weeks  that  might  be  required  to  alter  laboratory  apparatus  and  reduce  meas¬ 
urements  data.  The  program  is  worth  pursuing. 


PART  1 1 


ANALYTICAL  FORMULATION  FOR  EVALUATING 
INFRARED  SHIELD  CONFIGURATIONS 


SUtWARY 


An  analytical  approach  has  been  developed  for  predicting  directional 
radiant  intensity  from  shielded  infrared  sources.  In  this  approach, 
the  geometrical  and  radiative  characteristics  of  the  shield  and  source 
are  embodied  in  terns  of  the  integral  equations  for  radiative  transfer. 
The  solutions  of  these  integral  equations  served  as  the  basis  for  develop 
ing  a  prediction  formula  for  diffuse  and  for  specular  surfaces. 

The  prediction  formula  was  applied  to  diffuse  configurations,  both  with 
and  without  a  spherical  innerbody,  and  the  results  were  compared  with 
laboratory  measurements  of  a  similar  set  of  configurations.  Qualita¬ 
tive  agreement  was  found  to  exist  for  those  configurations  that  were 
directly  comparable. 


/ 
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INTRODUCTION 


This  part  of  the  report  is  concerned  with  the  development  of  an  analyt¬ 
ical  foundation  for  prediction  of  radiant  Intensity  patterns  from  shielded 
infrared  sources.  The  analytical  approach  roust  be  necessarily  founded 
upon  considerations  of  source  and  shield  geometry  and  of  the  radiative 
properties  of  the  surfaces.  It  is  the  object  of  this  report  to  incor¬ 
porate  these  basic  considerations  into  a  comprehensive  theoretical 
framework  from  which  the  shielding  effect  of  particular  configurations 
and  radiative  characteristics  can  be  derived. 

It  is  pertinent  to  mention  the  generality  with  which  the  infrared  shield¬ 
ing  problem  may  be  formulated.  From  the  geometry  aspect,  it  is  desirable 
to  have  a  theoretical  formulation  which  will  handle  a  wide  range  of  source 
and  shield  configurations.  On  the  other  hand,  among  the  radiative  prop¬ 
erties  of  the  source  and  shield,  the  general  formulation  should  account 
for  the  followingi 

(1)  Surface  distribution  of  radiant  emittance  or,  equivalently, 
the  surface  temperature  distribution. 

(2)  Directional  distribution  of  radiant  emissioni  i.e.,  the  re¬ 
lation,  for  an  elementary  unit  area,  of  the  radiant  intensity 
in  any  direction  to  that  in  the  normal  direction. 

(3)  Directional  distribution  of  reflected  radiant  energy  from 
an  elementary  unit  area,  accounting  for  the  specular  and 
diffuse  components  of  the  reflected  energy. 

(4)  Spectral  distribution  of  emitted  and  reflected  energy. 

(5)  Emissivity  and  reflectivity  of  the  source  and  shield 
materials. 

The  general  formulation  must  also  account  for  any  enclosed  attenuating 
media,  which,  while  they  are  not  properties  of  the  enclosing  surfaces, 
may  have  very  pronounced  effects  on  the  energy  transfer.  Furthermore, 
the  general  framework  should  handle  the  interreflection  of  radiation 
within  the  source-shield  configuration,  and  it  should  account  for  the 
presence  of  innerbodies.  This  ideal  generalization  is  of  more  than  aca¬ 
demic  value,  for  it  serves  as  a  standard  against  which  any  formulation 
must  be  compared.  It  will  be  shown  that  the  foregoing  generalizations 
of  the  problem  can  be  achieved  through  use  of  an  integral  equation  for 
radiant  energy  transfer.  This  implies  that  complicated  radiant  Intensity 
prediction  problems  can  be  formulated  and  solved,  in  principle  at  least, 
in  terms  of  and  through  the  use  of  integral  equations. 


But  there  arises  the  logical  question  of  to  what  extent  can  the  different 
aspects  of  the  problem  be  treated  simultaneously  In  general  terms.  This 
Is  a  question  of  great  practical  consequence  because,  for  a  given  source- 
shield  configuration,  simplifying  conditions  or  assumptions  are  usually 
required.  For  example,  the  condition,  or  assumption,  of  no  atmospheric 
attenuation  is  usually  taken  for  granted.  Another  frequently  imposed 
condition  is  that  of  gray,  diffusely  emitting  surfaces,  i.e.,  surfaces 
whose  emissivities  are  independent  of  wavelength  and  whose  directional 
distribution  of  emitted  energy  is  diffuse.  It  is  possible  to  meet  both 
these  conditions  to  a  fair  extent  over  spectral  bands  within  atmospheric 
windows,  and  both  conditions  are  assumed  in  this  report. 

The  next  section  presents  the  theoretical  considerations  used  in  formu¬ 
lating  the  shielding  problem.  The  integral  equation  for  radiant  energy 
transfer  is  derived  first  in  the  most  general  terms  and  then  is  particu¬ 
larized  to  apply  to  configurations  and  boundary  conditions  imposed  by  the 
experimental  work.  Three  basic  shield  configurations  are  considered; 
namely,  (1)  the  cylinder,  with  and  without  a  spherical  innerbody,  and 
(3)  the  concave,  truncated  double  cone.  Only  perfectly  diffuse  and 
specular  surfaces  are  considered.  In  the  third  section  are  presented 
graphical  comparisons  between  the  predicted  and  experimental  radiant 
intensity  patterns  for  diffuse  shields.  The  effects  of  individual  shield 
parameters  are  discussed,  and  suggestions  for  optimizing  shield  design 
are  made. 
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THEORETICAL  BASIS 


DERIVATION  OF  INTEGRAL  EQUATION 
Arbitrary  Enclosure 

Integral  equations  have  long  been  used  in  the  analysis  of  radiant  energy 
transfer  in  at  least  two  different  contexts,  that  of  illuminating  engi¬ 
neering  and  that  of  thermal  radiant  transfer.  Though  the  use  of  integral 
equations  in  these  fields  is  not  new,  their  usage  has  not  been  very  ex¬ 
tensive.  A  derivation  of  the  integral  equation  for  radiant  energy  trans¬ 
fer  is  given  here  to  illustrate  the  principle  and  to  emphasize  the  inherent 
generality.  The  derivation  follows  along  that  given  by  Reference  1. 


Figure  26.  Arbitrary  Radiating 
Enclosure. 


Consider  an  arbitrary  enclosure,  as  shown  in  Figure  26,  and  designate 
small  area  elements  at  two  arbitrary  points  by  dAi  and  dA.  The  positions 
of  the  two  area  elements  are  specified  in  terms  of  a  coordinate  system 
superimposed  on  the  surface  of  the  enclosure.  Suppose  that  element  dA^ 
is  held  fixed  at  coordinates  (x^,y^)  and  allow  element  dA  to  be  located 
at  the  variable  point  (x,y).  Attention  is  now  directed  to  the  fixed 
point  (xi,yi)  and  the  area  element  dA^,  at  which  a  radiant  energy  balance 
will  be  set  up.  The  total  radiant  energy  leaving  dAi  has  two  components. 
The  first  is  the  directly  emitted  energy,  which  in  the  absence  of  interre 
fleeted  radiation  is  called  the  emittance  Wo(xi,y^)  at  dA^.  The  second 
component  is  the  reflected  portion  of  the  irradiation  H(x^,yi)  of  dA]^ , 
which  is  given  by  ^(x^y^)  H(x1,yi),  where  />(xi,yi)  is  the  reflectivity 
of  the  enclosure  surface  at  (x]_,y]J.  Hence,  the  radiant  energy  balance 
at  dAj^  may  be  expressed  by 

*(xi,yi)  =  *o(xl*yl)  "r  l,yl)  H(xl»Yl)»  (3) 
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where  W(xj[,y^),  the  apparent  emittance  at  dA^,  and  H(xi,y^)  are  both  un¬ 
knowns.  The  term  H(xi,yi)  can  be  formally  eliminated  by  considering  the 
source  of  irradiation  of  dA^.  Thus,  consider  now  the  element  dA  at  any 
other  point  (  x,y  ).  If  the  total  radiant  energy  leaving  dA  is  W(x,y)» 
then  the  amount  reaching  dAi  is  equal  to 


dH(xi,yi)  =  W(x,y)  K(x,yi  xltyi)  dA,  (4) 


where  K(x,yj  x, ,Yi)  is  equal  to  the  irradiation  at  dA^  caused  by  unit 
(apparent)  emittance  from  unit  area  at  point  (x,y).  The  total  irradia¬ 
tion  of  dA^  is  obtained  by  integrating  the  contributions  from  all  the 
surface  elements  which  are  visible  from  dA^.  Thus  it  follows  that  Eq.  (3) 
can  be  restated  in  terms  of  a  single  unknown  function  as  follows: 


W(xi,y1)  =  W0(x1,y1)  +  p(x: ,yi )  K(x,yj  xltyi)  dA. 

A, visible 


(5) 


Since  the  unknown  function  W  appears  also  under  the  integral  sign,  Eq.  (5) 
is  an  integral  equation. 

Equation  (5),  which  may  be  called  the  integral  equation  for  radiant  energy 
transfer,  applies  to  an  arbitrary  enclosure  shape  and  to  arbitrary  surface 
radiative  properties.  That  this  is  true  is  seen  from  the  meaning  of  Eq.(4), 
for  it  contains  no  assumptions  at  all  on  the  geometry  nor  on  the  nature  of 
the  radiating  surface  element  dA.  It  merely  relates  the  total  energy 
W(x,y)  leaving  dA  in  all  directions  to  that  fraction  of  W(x,y)  which  im¬ 
pinges  upon  dAi,  all  other  circumstances  of  the  transfer  being  contained 
in  the  function  K.  Concisely,  the  function  K  contains  the  effects  of 
geometry,  of  the  directional,  and  spectral,  distributions  of  the  emitted 
and  reflected  radiation,  and  of  any  attenuating  media.  Thus  in  the  pres¬ 
ent  formulation  a  major  problem  revolves  around  the  specification  of  the 
function  K,  which  is  called  the  kernel  of  the  integral  equation.  A  more 
thorough  description  of  the  kernel  is  given  in  a  later  section. 


Cone  With  Spherical  Innerbody 

However,  before  considering  the  explicit  nature  of  the  kernel,  it  is  nec¬ 
essary  to  adapt  the  generalized  integral  equation  to  a  particular  geometry. 
For  this  purpose,  an  example  will  be  made  of  a  conical  tube  with  a  spheri¬ 
cal  innerbody  located  on  the  longitudinal  axis,  as  shown  in  Figure  44 
(Appendix  A).  Imagine  the  cone  and  sphere  to  be  divided  into  infinitesi¬ 
mal  rings,  or  bands,  transverse  to  the  longitudinal  axis  of  symmetry. 
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Now  consider  the  radiant  energy  balance  at  any  point  (x, fii)  on  the  coni¬ 
cal  band  element  specified  by  the  position  variable  x.  The  variable  x 
is  measured  along  the  axis  of  the  cone-sphere  system,  and  azimuthal  angles 
0  are  measured  from  the  upper  x-z  plane.  The  elementary  area  at  point 
(x,0^)  is  termed  a  "point”  element,  which  is  given  by 


dAj  =  r(x)  d  p  ^  dx  sec  0, 


where  r(x)  is  the  radius  of  the  cone  at  x,  fi.  Is  the  azimuthal  angle  of 
dA|,  and  0  is  the  cone  angle.  The  energy  balance  at  dAi  will  obviously 
involve  radiation  from  the  visible  portion  of  the  sphere  as  well  as  from 
the  remaining  unobscured  portion  of  the  cone  itself.  At  the  same  time, 
the  energy  balance  at  point  element  dA3  on  any  spherical  ring  will  involve 
radiation  from  the  unobscured  portion  of  the  cone  as  seen  from  dA3.  The 
point  element  dA3*  on  a  spherical  ring  specified  by  polar  angle  1  ,  is 
given  by 


dA3  =  R2  sin*  dr  d  03, 


where  R  is  the  radius  of  the  sphere  and  ^3  is  the  azimuthal  angle  of  dA3« 
Therefore,  the  total  energy  leaving  a  point  element  such  as  dAi  is  com¬ 
posed  of  two  parts,  the  first  being  the  intrinsic  thermal  emittance 
W0i(x,  0^)  at  the  point.  The  second  part  is  the  reflected  portion  of  the 
energy  received  from  elementary  ring  elements  of  the  cone  and  sphere  as 
seen  from  dAi.  Analogous  statements  can  be  made  for  the  spherical  point 
element . 

Before  it  is  possible  to  put  the  energy  balance  equation  into  useful  mathe¬ 
matical  form,  it  is  necessary  to  specify  the  obscuration  caused  by  the 
innerbody.  For  axial  symmetry,  the  obscuration  caused  by  a  spherical  inner- 
body  within  a  truncated  cone,  closed  at  one  end,  is  analyzed  in  Appendix  A. 
The  main  feature  of  this  analysis  is  that  the  obscuration  effect  is  accounted 
for  on  each  infinitesimal  ring  element. 

/ 

Because  of  the  axial  symmetry,  it  is  possible  to  express,  for  example,  the 
radiant  flux  from  a  partially  obscured  conical  ring  to  a  conical  point 
element  as  a  fraction  of  the  flux  from  the  complete  ring.  The  fraction 
varies  from  ring  to  ring  and  of  course  changes  for  any  given  ring  when  the 
observation  point  is  translated.  As  shown  in  Appendix  A,  these  facts  are 
embodied  in  a  function  T(x,f  ),  where  x  is  a  position  of  the  observation 
point  and  f  is  a  position  variable  also  along  the  axis  of  the  cone.  The 
function  T(x, f  )  is  bounded  by  zero  and  unity,  the  zero  value  correspond¬ 
ing  to  total  obscuration  and  the  unity  value  to  none.  On  the  other  hand, 
the  aximuthal  angular  limits  of  obscuration  for  each  spherical  ring  must 
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be  determined  for  each  and  every  observation  point  dA^  on  the  cone. 

Similar  limits  must  be  found  for  the  conical  rings  when  the  viewing 
point  is  taken  on  the  sphere.  Appendix  A  also  treats  the  obscuration 
of  the  end  disk,  which  is  introduced  later  below. 

With  the  obscuration  specified  according  to  Appendix  A,  attention  is 
turned  again  to  the  energy  balance  equations.  It  is  important  first  to 
point  out  that  sets  of  simultaneous  integral  equations  can  be  derived  on 
the  basis  of  the  energy  balance  for  the  ring  to  point  element  interchange. 
This  consideration  is  necessary  when  axial  symmetry  does  not  exist  or 
when  the  emittance  is  not  uniform  around  the  ring  element,  i.e.,  when  it 
is  a  function  of  the  azimuthal  angle  fi  .  For  the  cone-sphere  system  under 
consideration,  the  ring  emittance  is  assumed  to  be  independent  of  . 
Hence,  the  energy  balance  equation  at  point  element  dA^  with  a  conical 
band  at  f  and  a  spherical  ring  at  /'  can  be  written  as 


Wi(x)  =  W01(x)  +/>LW1(5)T(x,f)kn(x,f)da1  v(x,*0da3  (6a) 

where  da^  =  2W’£($)dE  is  the  area  of  the  conical  band  at  point  f  and 
da3  =  2«’R2sin7'd/'  Is  the  area  of  the  spherical  ring  at  The  analo¬ 

gous  equation  applying  at  the  spherical  point  element  dA3  is  given  by 


W3(*)  =  W03(*)  +^3  WlW  *13, v  dai’  (6b) 


where  in  both  equations  the  symbol  v  on  the  kernel  Indicates  that  only 
the  visible  portion  of  the  ring  is  considered.  The  total  energy  W]/  (x) 
at  dhi  follows  from  Eq.  (6a)  by  integrating  ^  over  tN»  length  of  the  cone 
and  f  over  the  visible  segment  of  the  sphere,  while  the  total  energy 
W^OO  at  dA3  follows  from  integrating  over  the  visible  portion  of  the 
cone.  Thus  the  following  equations  for  ring-to-point  interchange  result: 


i 


f* 2(x) 

*i(x)  =  w01(x)  */l|  /  *l(f  )f(x,f  Jknd*!  4  W^Qr0k31,v<J*3 


A/x) 


w3(*)  =  "03W 


(7a) 


l*f)  k13,vdal» 


fi  (r> 
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where  the  functional  notation  has  been  momentarily  suppressed  in  the  k 
functions.  (Subscript  1  denotes  the  cone  while  subscript  3  denotes  the 
sphere ). 

Up  to  this  point,  attention  has  been  concentrated  on  the  energy  balance 
at  only  point  elements.  With  the  symmetry  and  conditions  described 
above,  it  is  really  more  convenient  to  consider  the  receiving  area  ele¬ 
ment  to  be  a  ring.  The  energy  balance  equations  will  not  be  changed  in 
form,  since  the  energy  received  by  a  uniform  ring  element  is  just  a  con¬ 
stant  multiple  of  that  received  by  one  of  its  constituent  point  elements. 
Therefore,  on  the  basis  of  ring-to-ring  interchange  ,•  Eq.  (7a)  takes  the 
following  formt 


*l(x)  =  *0lM  + 


Wx(f  )  "1-  (x,f  )  ^n(x»? ) 


dr 


f*2(x) 

+1  w3(r)  K31tV  (x,y)  d7" 

^l(x) 

w3(r)  =  w03(t)  +  f  J  wx(f)  k13>v  (r, f)  dr, 

x(r) 


(7b) 


where,  for  convenience,  the  parameters  of  the  differential  ring  elements, 
except  the  differential,  have  been  combined  into  the  related  k  function. 

This  new  set  of  kernels  is  denoted  by  the  symbol  K.  Equation  (7b)  is 

the  desired  adaptation  of  the  generalized  form,  Eq.  (5),  with  two  integral 

equations  being  necessary  to  describe  the  conical  tube-sphere  configuration. 

This  special  case  has  been  considered  for  both  concreteness  and  simplicity 
'  in  order  to  illustrate  the  underlying  principles.  The  same  principles  and 
approach  are  applicable  to  all  configuration,  the  end  result  being  that  an 
integral  equation  is  required  for  each  separate  surface  of  the  enclosure. 

As  yet,  the  only  assumption  is  that  each  ring  element  must  have  uniform 
apparent  emittance  around  its  periphery,  but  nothing  has  been  assumed 
about  the  directional  distribution  of  the  energy  leaving  the  ring  element. 
For  reference  in  this  context,. and  for  completeness,  the  system  of  integral 
equations  which  obtain  when  a  radiating  disk  source  (denoted  by  subscript  2) 
is  used  to  close  one  end  of  the  cone  is  presented  here: 
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*o 


Wi(x)  =  W01(x)  Wi(f)Ttx,f)Ku(x,f)d?  +  JW2(/»)ttx,/)lC2i(x^)  4/ 

*o 

■OfjCx) 

«3(^')K31,v(x»^)  d*‘ 


^(x) 


(8a) 


W2(rO  =  W02^) 


|wx  (f  )T( r*,  0  Kx  2  (^* 


p^(r') 

)df  4  *3WK32#v<r,^dr 


^(r1) 


(8b) 


"3(*>  =  *03W+fi 


ri 


& 


f 

#o 


Wl(f)Kl3,v(r»^df  +  *2(f)K23,v(*'f)d/’ 


(8c) 


where,/*  and  r‘ are  radii  of  elementary  plane  annuli  on  the  disk  and  a 
is  the  radius  of  the  disk.  The  functions  7r(r,f)  and  T(*rf)  relating  to 
obscuration  of  the  end  disk  and  the  angular  limits,  oCl  and«2»  of  obscura¬ 
tion  pertaining  to  the  sphere  itself  are  given  in  Appendix  A. 


INTERPRETATION  OF  THE  KERNEL 


The  physical  content  of  thekernel  K  was  mentioned  briefly,  earlier. 
Attention  is  now  directed  to  the  inherent  significance  of  the  kernels  and 
to  establishing  a  more  thorough  physical  interpretation  of  them.  A  limited 
discussion  of  the  mathematical  nature  of  the  kernel  and  its  relation  to  the 
solution  of  the  integral  equation  is  also  included  in  this  section. 


Physical  Interpretation 

It  has  been  pointed  out  that  the  kernel*  of  the  integral  equation  repre¬ 
sents  a  major  portion  of  the  physical  effects  of  the  radiating  enclosure. 

*  It  should  also  be  pointed  out  that  the  kernel  as  used  here  has  been  called 
variously  the  angle  factor,  the  configuration  factor,  or  the  shape  factor. 
However,  these  names  often  relate  only  to  the  geoewtrical  nature  of  the 
problem. 
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Though  it  is  not  difficult  in  concept  to  include  atmospheric  attenuation 
in  a  physical  interpretation  of  the  kernel,  the  mathematical  details  be¬ 
come  exceedingly  complex.  The  inclusion  of  an  absorbing  medium  has  been 
treated  for  a  few  simple  geometries (Reference  2)  but  for  the  remainder 
of  this  report  the  effects  of  attenuating  media  will  not  be  considered. 
Therefore,  the  formulation  will  be  limited  to  "in  vacuo"  conditions, 
leaving  for  consideration  the  effects  of  the  radiative  surface  properties 
and  geometry. 

A  quite  significant  characteristic  of  a  radiating  surface  element  is  the 
distribution  in  space  of  the  energy  leaving  the  element.  The  directional 
distribution  must  be  considered  in  specifying  the  kernel  for  radiant  in¬ 
terchange  between  different  surface  elements.  For  an  elementary  plane 
area,  the  distribution  function  may  be  expressed  3S  the  ratio  (Reference  3) 
of  radiance  N(fl£)  in  the  direction  making  an  angle  d  with  the  normal  direc¬ 
tion  to  tbe  radiance  NQ  in  the  normal  direction,  or  as  the  ratio  (Ref¬ 
erence  4)  of  directional  emissivity  to  hemispherical  emissivity.  The 
first  ratio  is  just  ^ 


D(<*}  =  N(*)/N0. 


The  radiance  NQ  is  considered  to  be  a  mean  value  given  by  N  =  w/^f , 
where  W  is  the  radiant  emittance  of  the  surface  element.  TRe  distribu¬ 
tion  function  may  range  in  form  from  Lambert's  cosine  law  for  perfectly 
diffuse  surfaces  to  that  for  highly  directional  emission.  Lambert's  Law, 
namely, 


N(<*)  =  Nq  cosot  , 


is  frequently  assumed  for  the  emitted  radiation. 

For  a  cavity  surface  element,  the  directional  distribution  may  be  extended 
to  include  the  energy  reflected  from  the  surface  element,  since  the  com¬ 
ponent  of  reflected  energy  combines  with  that  directly  emitted  in  a  given 
direction  to  form  an  effective  total  intensity  in  that  direction.  The 
reflective  properties  of  the  surface  element  determine  in  part  the  direc¬ 
tional  distribution  of  the  energy  leaving  the  surface  element.  For  ex¬ 
ample,  a  diffusely  emitting  surface  may  be  a  specular  reflector;  or  it 
can  be  a  diffuse  reflector,  requiring  in  the  first  case  a  highly  polished 
surface  and  in  the  second  a  dull,  matte  finish.  In  the  second  case,  the 
resultant  distribution  would  remain  Lambertian,  whereas  any  other  combina¬ 
tion  would,  strictly  speaking,  result  in  a  non-Lambert ian  distribution. 
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Following  the  lines  of  Reference  3,  the  explicit  geometrical  character 
of  the  kernel  will  be  demonstrated.  For  this  purpose  consider  Figure  27, 
which  shows  the  geometrical  relation  between  two  elementary  areas  dA^  and 
dA2  separated  by  a  distance  r.  Let  o£^  be  the  angle  formed  by  the  normal 
direction  to  dA^  and  the  connecting  line  r  between  dA^  and  dA2,  and  let 
dw^  be  the  solid  angle  subtended  by  dA2»  Then,  that  energy  from  dAj  which 
impinges  on  dA2  is 


dE12  -  N^)  dwL  dA: 


(9a) 


whereas  the  total  energy  from  dAj  into  the  entire  hemisphere  is 


dE  =  dAL  =  TfNodA^  (9b) 

/ 

Hence,  by  definition,  the  kernel  K^2  is  given  by 

K12  =  dE^  =  >  “"I  .  (10a) 

dE  N0 


Since  the  solid  angle  dwi  is  given  by 


dwj 


dA2  coso£  2 


» 
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Equation  (10a)  can  be  expressed  as 


K12  =  Nl*!)  dA2  COSOt  2  _ 

n  N0  r2 


(10b) 


Equation  (10b)  represents 
tion  N(«0/No. 

Equation  (10b)  can  be  put 

dAx  K12 


the  kernel  for  an  unspecified  distribution  func- 

into  the  following  form: 

Dl«i)  coso(1  cos  #2  dA1  dA2  ,  . 

=  - - -  •  (Ha) 

ft  rl 


When  the  roles  of  dA^  and  dA2  are  reversed  in  the  foregoing  derivation, 
i.o.,  by  considering  dA2  to  be  the  emitting  surface  instead  of  the 
receiver,  an  equation  analogous  to  Eq»  (11a)  is  obtained,  i.e., 


dA2  K2 1 


P2(^2)  cosot  1  cos<T2  dki  dA2 

if  r2 


(lib) 


Equations  (11a)  and  (lib),  holding  for  any  directional  distribution  of  the 
corresponding  point  source  elements,  also  clearly  demonstrate  the  geometri 
cal  dependence  of  kernels  Kj2  and  K2^. 

It  is  clear  that  if  the  distribution  functions  (ot ^ )  and  ^(fl^)  are  equal 
then  equality  obtains  between  Eqs.  (11a)  and  (lib),  i.e., 

dAx  K12  =  dA2  K21.  (12) 


The  relation  expressed  by  Eq.  (12)  is  called  the  reciprocity  theorem. 
Reciprocity  of  the  radiant  interchange  between  two  area  elements  thus 
requires  identical  directional  distribution  at  each  element.  The  equality 
of  D(«fi)  and  0(0(2)  *s  assured  when  the  two  elements  are  Lambertian,  for 
then  D(*)  becomes  equal  to  unity  and,  hence,  independent  of  <L  .  The 
assumption  of  the  reciprocity  theorem  leads  to  a  very  useful  simplifica¬ 
tion,  since  it  effectively  reduces  by  one-half  the  number  of  required 
kernels  and  subsequent  computational  efforts. 

The  kernel  for  radiant  interchange  between  two  surface  elements  can  be 
derived  analytically  by  two  methods.  The  first  (References  3  and  5),  and 
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fundamental,  technique  is  integration  of  Eqs.  (11a)  or  (lib)  over  the 
appropriate  area  elements}  the  second  (References  3,  5,  6,  7,  8)  makes 
use  of  algebraic  operations  and  the  reciprocity  relation  to  derive  re¬ 
quired  kernels  from  those  already  existing.*  These  techniques  have  oeen 
applied  to  the  diffuse  shield  configuration  considered  in  this  report; 
the  results  are  compiled  in  Appendix  B. 


Mathematical  Classification 

Each  of  the  kernels  presented  in  Appendix  B  is  ?  function  of  the  posi¬ 
tions  of  two  differential  area  elements,  each  of  which  is  specified  by 
a  single  coordinate.  Always  one  point  is  fixed  while  the  other  varies 
over  an  open  interval.  The  functional  behavior  of  the  kernel  K(x,y)  is 
of  major  importance  in  the  formal  theory  of  integral  equations.  Without 
pretense  of  more  than  just  cursory  mention,  some  of  the  mathematical 
characteristics  of  the  kernel  are  listed  here.  The  kernel  K(x,y)  is 
required  to  be  real  but  not  Identically  zero  on  the  interval  of  interest. 
Theory  permits  that  both  K(x,y)  and  the  yet  unknown  function  W(x)  may  be 
discontinuous  functions,  provided  they  are  of  integraole  square.  Two  im¬ 
portant  classifications  exist  which  depend  on  the  functional  form  of  the 
kernel.  First,  the  kernel  is  of  the  separable  type  when  K(x,y)=h(x)g(y) , 
or  a  sum  of  such  products,  and  the  second  type  is  the  symmetric  kernel, 
i.e.,  when  K(x,y) =  K(y,x).  Their  consequence  in  relation  to  the  solution 
of  the  integral  equation  is  pointed  out  in  the  next  section.  For  the 
symmetric  kernel,  a  rather  complete  description  of  mathematical  theorems 
and  details  and  their  relation  to  interref lecting  enclosures  has  been 
presented  in  Reference  9. 


SOLUTION  OF  THE  INTEGRAL-EQUATION  SYSTEM 

The  approach  for  solving  the  integral-equation  system  is  determined  by  the 
nature  of  the  kernel  involved.  In  turn,  the  kernel  is  directly  dependent 
on  the  symmetry  of  the  shield  configuration.  An  approximate  analytical 
solution  is  possible  for  the  cylindrical  shield,  under  certain  conditions, 
whjlle  a  direct  numerical  solution  must  be  obtained  for  the  other  configura¬ 
tions  considered  in  this  report.  The  conditions  under  which  (approximate) 
analytical  solutions  obtain  are  presented  below,  followed  by  a  description 
of  the  numerical  method  of  solution. 


*  The  references  cited  here  on  this  subject  are  intended  to  be  repre¬ 
sentative  rather  than  exhaustive.  Reference  6  gives  a  review  of  defi¬ 
nitions  and  methods  for  deriving  kernels  (or  angle  factors). 
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Analytical  and  Approximate  Analytical  Methods 


If  the  shield  itself  is  sufficiently  symmetrical,  then  the  corresponding 
kernel  is  synmetric.  For  a  surface  of  revolution,  at  least  three  axes 
of  symmetry  are  required  to  assure  a  symmetric  kernel.  The  Hilbert -Schmidt 
method  (see  Reference  10,  for  example)  applies  to  the  case  of  a  symmetric 
kernel,  allowing  in  principle  closed-form  solutions.  But  even  symmetric 
kernels  lead  to  integrals. which  are  not  readily  integrated.  Such  cases 
are  the  rule,  and  it  becomes  necessary  to  replace  the  true  kernel  with 
a  simpler  analytic  approximation.  The  usual  stipulation  on  the  form  of 
the  approximate  kernel  is  that  it  be  separable.  The  condition  of  separa¬ 
bility  results  in  simpler  integrals  which  can  then  be  integrated  to  give 
ar.  approximate  solution. 

It  has  been  found  (Reference  11)  that  rather  simple  approximate  kernels 
give  satisfactory  results  for  a  cylindrical  enclosure.  The  approximate 
kernels  are  in  the  form  of  a  two-  or  three-term  exponential  function. 

Even  a  one-term  exponential  has  been  found  by  Moon  (Reference  12)  to  be 
satisfactory  for  illumination  calculations.  Moon's  solution  is  based  on 
the  assumption  that  the  initial  and  final  emittances  are  uniform  ever  the 
end  disk.  His  solution  was  applied  to  the  cylinder  in  the  present  case, 
and  the  result  is  compared  with  a  numerical  solution  in  the  next  section. 

The  foregoing  discussion  points  to  the  fact  that  even  approximate  analytic 
solutions  seem  possible  for  only  highly  symmetrical  shield  configurations. 
Such  analytic  solutions,  however,  possess  the  inherent  advantage  and  gen¬ 
erality  of  closed  form,  in  which  the  effects  of  the  separate  parameters 
of  the  shield  can  be  examined  with  ease;  but  the  results  are  approximate 
because  the  true  kernel  has  been  replaced  by  an  approximation.  For  a 
wider  class  of  less  symmetrical  shields,  the  approach  rrwst  be  along  the 
lines  of  numerical  solutions. 


Numerical  Solution 


Direct  numerical  solutions  have  one  distinct  disadvantage  in  that  they 
, have  no  generality  of  form  and  must  be  worked  out  again  for  each  combi¬ 
nation  of  parameters.  In  effect,  however,  this  disadvantage  can  be  made 
almost  negligible  through  the  use  of  high-speed  electronic  computers.  In 
addition,  numerical  solutions  can  be  achieved  for  a  precise  and  complete 
formulation  of  the  physical  processes  within  the  shield  without  resort  to 
approximate  kernels.  This  fact  is  especially  important  when  a  system  of 
integral  equations,  such  as  Eq.  (8),  is  required,  for  then  approximate  an¬ 
alytical  solutions  would  be  difficult,  if  possible  at  all.  This  state  of 
affairs  may  be  summarized  as  follows:  Approximate  numerical  solutions  are 
more  readily  obtained  for  a  complete,  analytical  formulation  of  the  shield 
problem  than  are  analytical  solutions  for  an  approximate  formulation. 
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The  numerical  solution  technique  that  was  used  in  the  present  work  is 
now  described.  Similar  and  other  methods  of  solution  have  been  described 
In  the  literature  in  relation  to  other  problems  involving  diffuse  (Ref¬ 
erences  7  and  8)  and  specular  (References  13  and  14)  enclosures.  The 
method  of  solution  used  herein  is  an  iteration  technique,  which  applies 
to  either  diffuse  or  specular  shields.  Then,  considering  Eq.  (8),  the 
solution  procedure  is  outlined  as  follows:  First,  the  physical  problem  was 
specified,  for  any  given  geometrical  configuration,  by  assigning  appro¬ 
priate  values  to  the  radiative  parameters,  i.e.,  to  the  reflectivities 
and  the  initial  emlttances.  Next,  initial  trial  solutions  —  designated 
by  W^il)  and  W^f^l)  —  were  guessed  for  the  unknown  functions 
Wj(f )  and  Wt( Y’)  ,  respectively,  in  Eq.  (8b).  Ultimately,  each  initial 
trial  function  took  the  form  of  a  finite  sequence  of  numbers  dependent  on 
the  corresponding  position  variable  f  or  ;  and  each  was  introduced  as 
such  into  the  corresponding  integral  of  Eq.  (8b).  A  trial  distribution 
was  then  obtained  for  W2(r')by  numerically  evaluating  the  right-hand  side 
of  Eq.  (8b)  for  each  value  of  r',  the  end  disk  being  divided  into  a  finite 
number  of  evenly  spaced  radial  increments.  Let  W2(r';l)be  the  solution 
thus  obtained.  Then  the  new  function  W2(r';l)and  initial  guessed  function 
WiGfjl)  were  applied  in  Eq.  (8c),  and  in  a  similar  fash'.on  a  new  distri¬ 
bution  W3(3T|2)  was  generated  for  the  sphere.  Both  W3(fc'i;2)  and  K^r’jl) 
were  then  substituted  in  Eq.  (8a),  and  the  right-hand  side  was  numerically 
evaluated  for  each  value  of  x,  giving  a  second  trial  distribution  W^(x|2) 
applying  to  the  conical  walls. 

Starting  again  with  Wi(Jp2)  and  W3()',|2)  in  Eq.  (8b),  the  entire  process 
was  repeated. 

The  foregoing  procedure  was  iterated  until  convergence  was  achieved  on  all 
three  functions.  The  convergence  criterion  was  set  as  a  specified  differ¬ 
ence  between  the  values  of  the  function  on  the  nth  and  (n-l)st  iteration. 
In  this  fashion,  any  desired  degree  of  accuracy  may  be  acquired. 

Obviously,  when  the  innerbody  is  absent,  obscuration  effects  are  absent, 
and  the  computational  effort  is  considerably  reduced.  When  consideration 
of  the  innerbody  is  eliminated,  Eq.  (8)  reduces  to  two  integral  equations, 
which  are  also  solved  by  the  iteration  method. 


PREDICTION  OF  DIRECTIONAL  RADIANT  INTENSITY  PATTERNS 


Up  to  this  point,  attention  has  been  concentrated  on  the  physical  occur¬ 
rences  within  the  radiating  enclosure.  This  has  been  necessary  to  provide 
a  general  framework  within  which  the  prediction  of  radiant  intensity  direc¬ 
tional  patterns  in  space  can  be  formulated  for  a  wide  range  of  conditions. 
It  is  possible  to  develop  simpler  prediction  schemes  without  the  compli¬ 
cations  of  integral  equations  if,  for  example,  interref lected  radiation 
Is  assumed  to  be  unimportant;  or  if  the  distribution  of  emittance  over 


all  surfaces  of  the  enclosure  remains  uniform  in  the  presence  of  diffuse 
reflection  (Reference  15).  A  simpler  formulation  is  desirable  if  the 
actual  physical  problem  permits  it.  However,  for  many  practical  cases, 
such  simplifying  assumptions  are  not  justified,  requiring  a  more  thorough, 
though  more  complicated,  description. 

Attention  is  now  turned  to  the  role  played  by  the  solutions  of  the  inte¬ 
gral  equation  system  in  formulating  the  prediction  problem.  As  will  be 
shown,  the  approach  to  the  prediction  formulation  depends  heavily  on  the 
nature  of  the  shield  surface;  and  three  distinctions  are  readily  made 
according  to  whether  the  surface  is  diffuse,  specular,  or  intermediately 
non-Lambert ian  between  these  two  extremes. 


Role  of  Integral-Equation  System  Solution 

The  solutions  of  the  integral-equation  system  pertaining  to  a  shielded 
source  are  not  themselves  solutions  to  the  directional  intensity  pre¬ 
diction  problem.  They  are  necessary  intermediate  steps  in  a  general 
treatment,  however.  It  will  be  iterated  that,  physically,  the  solutions 
of  the  integral-equation  system  represent  the  distribution  of  apparent 
emittance  over  the  surfaces  of  the  shield-source  enclosure.  The  emit- 
tance  distribution  function  thus  obtained  contains  the  increase  caused 
by  interre fleeted  energy  above  the  intrinsic  emittance  of  a  surface  ele¬ 
ment  at  the  existing  temperature.  Also,  the  apparent  emittance  distri¬ 
bution  can  contain  the  effects  of  temperature  gradients.  For  shields 
which  are  surfaces  of  revolution,  the  integral  equation  approach  readily 
handles  temperature  gradients  along  the  longitudinal  axis  of  symmetry, 
i.e.,  the  temperature  varying  from  ring  to  ring  but  remaining  constant 
around  the  circumference  of  each  ring  element.  An  asymmetrical  tempera¬ 
ture  distribution  could  be  handled,  but  with  more  difficulty. 

The  apparent  emittance  distribution  is  therefore  quite  important,  for  it 
specifies  the  total  amount  of  radiant  energy  leaving  each  surface  element 
of  the  shield-source  enclosure.  The  apparent  emittance  of  an  elementary 
unit  area  gives  the  total  energy  transferred  therefrom  into  a  hemisphere 
and  thus  does  not  regard  the  directional  distribution  of  the  energy,  h'ow- 
'ever,  it  should  be  noted  again  that  if  the  total  energy  leaving  the  unit 
area  remains  diffuse,  then  its  directional  distribution  will  be  given  by 
Lambert's  cosine  law.  Such  is  not  the  case  for  specular  and  other  non- 
Lambert  ian  surfaces.  Whatever  the  case,  the  apparent  emittance  distribu¬ 
tion  serves  as  the  starting  point  for  the  prediction  formulations,  which 
are  now  considered. 


Formulations  for  the  Diffuse  Shield-Source  Enclosure 


Strictly  speaking, the  geometry  and  solid  angle  of  the  detector  optical 
system  relative  to  the  shield  aperture  are  important  factors  to  be 
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considered  in  deriving  an  intensity  prediction  theory.  This  is  especially 
pertinent  for  detector-source  separation  distances  such  that  the  field 
of  view  of  the  detector  optical  system  is  less  than  the  aperture  of  the 
shield  (Reference  16).  The  following  formulation  assumes  that  the  shield 
aperture  is  less  than  the  field  of  view  of  the  detector  and  that  the  col¬ 
lection  optics  receives  an  essentially  non-divergent  beam  of  radiation. 

It  is  well  known  that  the  radiant  intensity  JQ  in  the  direction  normal  to 
a  uniform  diffuse  source  is  simply  related  to  the  emittance  by 


J0  =  (1/4T)  WA, 


where  A  is  the  area  of  the  source.  If  the  emittance  W  were  a  function  cf 
position  on  the  source,  then  an  integration  of  W  over  the  source  would  be 
required  for  the  radiant  intensity,  i.e., 


J 


o 


W(x,y)dA. 


The  directional  distribution  of  radiant  intensity  J(0)  for  this  source 
would  follow  by  integrating  the  non-uniform  emittance  over  the  projected 
area  of  the  source  as  seen  from  an  aspect  angle  0,  treating  the  projected 
area  as  a  non-uniform  plane  source.  An  application  of  this  idea  was  made 
on  the  diffuse  cylinder  without  an  innerbody  under  the  assumption  that  the 
apparent  emittance  over  the  source  disk  remained  uniform.  However,  this 
method  was  difficult  to  apply  to  other  shield  configurations;  its  applica¬ 
tion  to  the  cylinder  will  be  described  later.  But  first  a  formulation 
applicable  to  all  the  shield  geometries  will  be  described. 

Consider  a  diffuse  conical  shield  with  diffuse  spherical  innerbody  and 
a  diffuse  circular  source  which  forms  the  closed  end  of  the  shield-source 
enclosure.  The  distribution  of  apparent  emittance  over  each  of  the  three 
surfaces  is  given  by  solving  Eq.  (3).  Generally,  the  apparent  emittance 
distribution  on  each  surface  will  be  non-uniform.  The  confronting  problem 
is  to  integrate  these  three  non-uniform  distributions  over  the  visible  por¬ 
tions  of  the  corresponding  surfaces  as  seen  from  different  aspect  angles. 
(See  Figure  28.) 

An  integration  in  closed  form  is  impossible,  if  for  no  other  reason  than 
the  fact  that  the  apparent  emittance  distribution  functions  are  not  in 
closed  form.  Hence,  a  finite  summation  process  must  be  used.  To  facili¬ 
tate  the  summation,  the  projected  area  defined  by  the  shield  aperture  at 
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the  desired  aspect  angle  0  is  considered  to  be  c  non-uniform  plane  source. 
This  plane  source  consists  of  the  projections  of  the  visible  portion  of 
the  interior  shield  wall)  the  innerbody,  and  the  source  disk.  As  shown 
in  Figure  28,  each  visible  surface  is  divided  into  a  number  of  small  area 
elements  which  are  bounded  by  elliptical  arcs.  This  latter  condition 
arises  from  the  fact  that  on  each  surface  the  apparent  emittance  is  cal¬ 
culated  at  the  centers  of  a  finite  number  of  circular  ring  elements j  these 
then  appear  as  elliptical  elements,  or  fractions  thereof,  in  the  projec¬ 
tion.  Each  elliptical  area  element  thus  can  be  weighted  by  the  apparent 
emittance  calculated  for  the  center  of  the  parent  circular  element.  The 
weighted  area  elements  are  then  summed  for  all  the  surfaces  visible  at 
the  given  aspect  angle  0  to  give  the  radiant  intensity  J(0)  in  that 
direction. 

A  mathematical  statement  of  the  foregoing  procedure  takes  the  following 
form: 


n  m 

J(0)  =  Ei  Alj(0,xj)W(xj)  +  C  i  ^k^’VW  * 


J2  *31(0.4>*3(4)  da) 


In  Eq.  (13),  the  A  A  terms  represent  the  elliptical  area  elements  de¬ 
scribed  above  on  the  projection  plane  at  aspect  angle  0,  where 


A.  (0,x  )  =  area  element  relating  to  the  circular  ring 
^  ^  element  centered  about  Xj  on  the  shield  wall. 

A2J((0,r  )  =  area  element  relating  to  the  annular  ring 

*  centered  about  radius  r ^  on  the  source  disk. 


A3i(0,Vi)  =  area  element  relating  to  the  spherical  ring 
centered  about  the  polar  angle  on  the 
spherical  innerbody. 

The  W  terms  are  the  corresponding  apparent  emittance  values.  The  sum 
involving  either  the  source  disk  or  the  innerbody  disappears  when  either 
is  completely  obscured  from  view.  The  analysis  establishing  the  three 
&  A  terms  as  a  function  of  angle  0  is  a  straightforward  application  of 
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Figure  9*  Interior  Surface  Elements  of  Conical  Enclosure 
with  Spherical  Innerbody. 


analytic  geometry;  these  derivation^  are  omitted  from  the  report.  Equa¬ 
tion  (13)  was  programmed  in  the  functional  forms  of  the  A  A  terms  for 
machine  computation  of  J(0). 

As  previously  mentioned,  a  second  prediction  method  can  be  formulated 
for  the  diffuse  shield-source  enclosure  when  the  apparent  emittance  dis¬ 
tribution  functions  are  known  in  closed  form.  Such  is  the  case  for  the 
cylinder.  The  approximate  analytical  solution  obtained  by  Moon  (Ref¬ 
erence  12)  for  the  diffuse  cylinder  is  of  the  following  form: 


W^(x)  =  Bcoshkx  +  Csinhkx 
W2 ( r  * )  —  constant 

For  the  experimental  conditions  considered  herein, 
duces  to  a  single  exponential  term  given  by 

W:(x)  =  Be"kx 

w2(rl)~  constant 


(14) 


Moon's  solution  re- 


(15) 


The  details  of  deriving  the  prediction  formulas  making  use  of  Eq.  (lb)  are 
given  in  Appendix  C.  As  shown  there,  two  prediction  formulas  are  required, 
depending  on  whether  the  source  disk  is,  or  is  not,  seen  at  a  given  aspect 
angle.  The6e  two  formulas  are  given  by  Eqs.  (51)  and  (52),  Appendix  C. 


Formulation  for  Specular  Shield 

Because  of  the  mirror-like  surface  elements,  many  aspects  of  the  specular 
shield  are  brought  within  the  realm  of  geometrical  optics.  The  predic¬ 
tion  formulation  becomes  essentially  a  ray-tracing  problem  of  varying 
difficulty,  depending  on  the  geometry  of  the  shield.  The  prediction  for¬ 
mulation  for  specular  shields  will  be  limited  herein  to  the  cylinder  closed 
at  one  end  by  a  diffuse  disk  source.  In  principle,  this  formulation  can 
be  extended  to  the  cone  and  other  geometries. 

It  is  necessary  first  to  determine  the  kernels  for  transfer  of  specular 
radiation.  The  kernel  for  specular  interchange  between  elementary  rings 
on  the  cylinder  is  given  by  Eq.  (8)  of  Reference  13.  The  kernel  for  specu¬ 
lar  interchange  from  an  elementary  annulus  to  an  elementary  ring  element 
on  the  cylinder  can  be  derived  by  the  method  of  Reference  13,  and  the  result 
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is  given  in  Appendix  B.  Each  of  these  specular  kernels  contains  a  con¬ 
vergent  infinite  series  which  represents  the  infinite  number  of  reflec¬ 
tions  possible  between  two  area  elements.  Next,  the  two  sinwltaneous 
integral  equations  follow  from  Eq.  (8)  by  eliminating  consideration  of 
the  innerbody.  Introducing  the  specular  kernels  and  initial  trial  solu¬ 
tion  into  the  integral  equations  and  applying  the  iteration  process  yields 
ultimately  the  composite  energy/unit  area  leaving  from  each  ring  element 
into  all  directions. 


It  is  now  possible  to  develop  a  prediction  technique  by  considering  the 
radiative  transfer  from  each  ring  element  to  the  openend  of  the  cylinder. 
As  an  example,  consider  a  ring  at  x  on  the  cylinder  in  Figure  29.  Radiant 
energy  from  the  ring  can  pass  through  the  open  end  either  directly  or  in¬ 
directly  by  reflection.  The  open  end  can  be  regarded  as  a  disk  which  is 
itself  divided  into  annular  zones  or  rings.  Each  annular  zone  in  turn  can 
be  thought  of  as  defining  a  small  solid  angle  into  which  the  entire  cylin¬ 
drical  ring  transfers  directly  a  fraction  of  its  total  energy.  From 
Figure  29,  the  direction  0(x)  of  this  small  solid  angle  with  respect  to 
the  axis  of  the  cylinder  is  given  by 


0(x)  =  arc  tan 


whereas  that  for  the  annular  ring  on  the  source  disk  is  given  by 


0(r  )  =  arc  tan 
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All  that  remains  to  determine  the  directly  transmitted  radiation  is  to 
specify  the  frac'on  of  the  total  energy  at  each  band  transferred  through 
each  zone  of  the  open  end.  It  turns  out  that  the  kernels  for  ring-to-ring 
interchange  multiplied  by  the  width  of  the  appropriate  annulus  or  ring 
are  the  required  fractional  values,  or  exchange  factors. 

It  is  possible  to  derive  an  exchange  factor  for  radiation  reaching  the 
open  end  via  n  reflections  (Reference  13).  The  angle  of  incidence  at 
the  first  reflection,  and  also  its  complement  0n ,  is  preserved  at  all 
n  reflections  because  of  the  cylindrical  symmetry!  the  angle  0n  is  the 
aspect  angle  of  the  n-tlmes  reflected  ray  with  respect  to  the  axis  of 
the  cylinder. 

From  the  ring  element  at  x,  for  example,  both  the  energy  transferred 
through  the  open  end  and  the  aspect  angle  0n  can  be  computed  for  once- 
reflected  rays,  for  twice-reflected  rays,  for  Nth  and  so  on  up  to  Mth- 
reflected  rays.  The  number  M  is  specified  by  the  criterion  that  the 
(M  +  1 )th-reflected  ray  transfers  negligible  energy  through  the  open  end. 
Similar  calculations  can  be  made  for  all  ring  elements  at  which  the  appar¬ 
ent  emittance  has  been  previously  calculated!  the  aspect  angles  0n(O  and 
0n(r')  and  the  corresponding  emergent  radiant  intensities  are  computed  and 
compiled  for  all  ring  elements.  From  the  total  collection  of  these  com¬ 
putations,  the  angles  0n(x)  and  0n(r’)*re  sorted  according  to  whether  they 
fall  within  a  neighborhood  A0  of  a  specified  aspect  angle  0.  The  radiant 
intensities  corresponding  to  the  emergent  directions  0n  falling  within 
0  +  A0  are  summed  to  give  the  radiant  intensity  of  the  reflected  energy 
in  the  direction  0. 

A  similar  sorting  technique  can  be  applied  to  the  directly  transmitted 
energy.  The  prediction  formulation  for  the  specular  cylinder  is  com¬ 
pleted  by  now  adding  the  direct  and  the  reflected  components  in  the  direc¬ 
tion  0. 
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COMPARISON  OF  PREDICTED  DIRECTIONAL  INTENSITIES 


The  results  of  calculations  based  on  the  prediction  formulations  devel¬ 
oped  in  the  preceding  section  are  now  presented  and  compared.  The  com¬ 
putational  results  are  presented  in  graphical  form  in  Figures  30  -43. 
Consideration  of  the  basic  shield  configuration  leads  to  four  logical 
groupings  which  are  keyed  to  tables  in  the  text.  The  corresponding 
table  specifies  all  the  other  shield  parameters  for  each  curve  of  the 
basic  group,  while  the  parameter  being  marled  is  listed  on  the  graph. 
Curves  not  in  these  groups  are  specified  within  the  figure.  A  direct 
comparison  with  laboratory  measurement  is  possible  for  a  set  of  param¬ 
eters  for  each  basic  shield  configuration. 

Description  of  the  experimental  work  is  given  in  the  interim  report  re¬ 
lating  .o  this  study  (Reference  18).  The  emlttance  of  the  unshielded 
source  based  on  the  experimental  work  is  .30  watts/cm2 j  this  value  was 
assigned  to  the  parameter  Wq2  in  all  calculations.  The  source  reflec- 
tivity/2  was  n°t  varied  either,  being  fixed  at  .30.  Since  the  shield 
models  were  maintained  at  room  temperature,  a  direct  measurement  of  the 
wall  emlttance  Wq^  was  not  possible  with  the  experimental  arrangement 
described  in  Reference  18.  Therefore,  an  extrapolated  value  of 
10~4  watts/cm2  was  taken  as  the  experimental  value  of  WqI  for  all  con¬ 
figurations.  The  reflectivity  of  the  experimental  diffuse  shiald  models 
was  .02.  The  parameters  WQ1  and  A,  for  given  values,  were  assumed  to 
be  uniform  over  the  shield  surface.  Attention  is  now  turned  to  comparison 
of  directional  intensity  patterns  and  the  effects  of  shield  paraiseters. 


EFFECTS  OF  SHIELD  PARAMETERS 
Diffuse  Cylinder 

Consider  first  the  diffuse  cylinder.  Directional  radiant  intensity  J(0) 
patterns  are  plotted  for  different  combination  of  parameters  in  Figures 
30,  31  and  32.  Figure  30  shows  the  dependence  of  J(0)  on  the  reflec¬ 
tivity  of  the  cylinder,  while  Figure  31  shows  the  dependence  on  the 
emlttance  of  the  cylinder.  Figure  32  gives  the  directional  patterns  for 
variation  in  the  cylinder  length,  and  also  compares  the  experimental  curve 
with  the  predicted  pattern  for  the  laboratory  conditions.  The  differences 
between  these  two  particular  curves  will  be  discussed  below  along  the  simi¬ 
lar  curves  for  the  other  configurations.  It  is  noted  from  Figure  30  that 
the  pattern  increases  with  increasing  reflectivity,  especially  for  the 
larger  values  of  0.  Since  the  wall  emlttance  is  held  fixed,  this  effect 
must  be  due  entirely  to  Interreflected  radiation  when  the  source  is  not 
visible,  which  occurs  for  this  particular  cylinder  when  0  is  greater  than 
27.11°.  A  distinct  bend  in  curves  2  and  3  starts  approximately  at  this 
angle,  indicating  clearly  the  increase  due  to  interreflection  over  the 


61 


0*001 


62 


Figure  30.  Directional  Intensity  of  Figure  31.  Directional  Intensity  of 

Diffuse  Cylinders  (Sep  Diffuse  Cylinders  (See 

Table  ‘O.  Table  b). 


Figure  32.  Directional  Intensity  of  Figure  33.  Directional  Intensity  of 

Diffuse  Cylinders  (See  Diffuse  Small  Cones  (See 

Table  b).  Table  6). 


directly  emitted  energy.  This  increase  can  be  taken  approximately  as  the 
difference  from  curve  1,  which  in  the  normal  direction  is  about  a  35%  in¬ 
crease  for  curve  4.  The  curves  of  Figure  31  show  the  expected  result 
that  as  the  cylindrical  cavity  approaches  the  isothermal  condition  the 
directional  intensity  pattern  tends  toward  a  cosine  distribution.  Figure  32 
clearly  shows  that  decreasing  the  cylinder  length  decreases  its  shielding 
effect. 


TABLE  5 

DIFFUSE  CYLINDER  PARAMETERS 


Fig.  No. 

Curve 

No. 

Length 

(in.) 

Radius 
(in. ) 

Cyl.  Ref 1 . 

fl 

WatU/cm* 

30 

1 

7.8125 

2 

.02 

.1  x  10“3 

- 

2 

7.8125 

2 

.30 

.1  x  J0“ 3 

- 

3 

7.8125 

2 

.60 

.1  x  10"3 

- 

4 

7.8125 

2 

.90 

.1  x  10r3 

31 

1 

7.8125 

2 

.02 

.1  x  10“3 

- 

2 

7.8125 

2 

.02 

.15 

- 

3 

7.8125 

2 

.02 

.30 

32 

1 

8.0 

2 

.02 

.1  x  10"3 

- 

2 

4.0 

2 

.02 

.1  x  10-3 

- 

3 

2.0 

2 

.02 

.1  x  10‘ 3 

* 

exp.* 

7.8125 

2 

.02 

.1  x  10'3 

*  Represents  experimentally  determined  curve. 

Diffuse  Small  Cone 


By  the  term  "small  cone"  is  meant  the  truncated  cone  with  its  wider  base 
forming  the  source  area.  The  curves  pertaining  to  this  configuration  are 
presented  by  Figures  33,  34,  and  35.  The  experimental  curve  for  the 
laboratory  small  cone  is  compared  with  the  corresponding  predicted  pat¬ 
tern  in  Figure  33.  A  comparison  between  curves  1  and  2  of  Figure  33 
shows  the  effect  of  changing  the  cone  angle  for  the  same  axial  length. 

The  directional  patterns  here  demonstrate  rather  abrupt  decreases  in 
radiant  intensity.  This  sharp  change  is  noted  for  the  experimental 
curve  also.  These  curves  follow  very  closely  the  cosine  law  just  short 
of  the  decrease.  Considering  curve  2,  Figure  34  (6  =  12.67°)  tad  curve  2, 
Figure  33  (9  =  11.31°)  shows  for  approximately  the  same  cone  angle  but 
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half  the  axial  length  that  the  cutoff  is  not  nearly  as  sharp.  This  is  to 
be  expected  since  the  open  radius  of  the  small  cone  increases  by  decreas¬ 
ing  the  axial  length  for  a  given  cone  angle.  This  effect  is  amplified  in 
Figure  35,  which  compares  several  cone  angles  for  an  axial  length  of  two 
inches.  The  isothermal  case,  represented  by  curve  4,  Figure  35,  results 
in  very  nearly  a  cosine  distribution,  indicating  little  shielding  effect. 


TABLE  6 

DIFFUSE  SMALL  CONE  PARAMETERS 


Fig. 

No. 

Curve 

No. 

Axial 

Length 

(in.) 

Base 
Rad. 
(in. ) 

Open 
Rad . 
(in.  ) 

Cone 
Ang. 
(deg. ) 

Cone 
Ref  1 . 

A 

W0i 

(watts/ca.^ ) 

33 

1 

8 

2.5 

1.375 

.8° 

.02 

.1 

x 

10“3 

- 

exp.* 

8 

2.5 

1.375 

.8° 

.02 

.1 

X 

10-3 

- 

2 

8 

2.5 

.90 

11.31° 

.02 

.1 

X 

10“3 

34 

1 

4 

2.5 

2.10 

5.71° 

.02 

.1 

X 

10*3 

- 

2 

4 

2.5 

1.60 

12.67° 

.02 

.1 

X 

10“3 

35 

1 

2 

2.5 

2.40 

2.86° 

.02 

.1 

X 

10*3 

_ 

2 

2 

2.5 

2.28 

6.41° 

.02 

.1 

X 

10-3 

- 

3 

2 

2.5 

1.58 

24.8  0 

.30 

.3 

X 

10~3 

v  Represents  experimentally  determined  curve. 


Diffuse  Large  Cone 

The  term  "large  cone”  means  the  truncated  cone  with  the  smaller  end  at  the 
source,  as  shown  by  Figure  44.  The  effects  of  varying  the  shield  param¬ 
eters  of  this  configuration  are  represented  by  Figures  3<-  through  40.  In 
Figure  36  are  plotted  the  directional  intensities  as  a  function  of  reflec¬ 
tivity  for  a  cone  angle  of  39. TO0.  These  curves  exhibit  a  sharp  cutoff 
similar  to  the  small  cone,  though  occurrin:  at  larger  aspect  angle  '.  The 
distribution  up  to  the  cutoff  angle  remains  almost  a  cosine  d istriojt ion , 
as  was  observed  experimentally.  Here,  as  in  the  cylinder,  the  effects  of 
interref lected  radiation  show  up  markedly  past  the  cutoff  at  which  the 
source  is  no  longer  visible.  Figure  37  ^hows  two  curves  applying  to  the 
same  cone  angle  of  68°,  approximately,  tut  for  lengths  of  two  and  four 
inches.  The  curves  coincide  almost  over  the  entire  quadrant  and  fall  in 
the  proximity  of  the  cosine  distribution.  Figure  38  demonstrates  the 
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(root  I - 1 - 1 - £=5— - -JOMI 


Figure  3B,  Directional  Intensities  of  Table  39.  Directional  Intensities  of 

Diffuse  Large  Cones  (See  Diffuse  Large  Cones  (See 

Table  7).  Table  7). 


variation  caused  by  changing  the  cone  emittance  #qi*  Figure  39  shows 
the  effect  of  changing  the  cone  angle.  It  seems  from  these  curves  that 
the  most  effective  shield  is  obtained  with  small  cone  angles,  approach¬ 
ing  the  cylinder.  The  predicted  curve  shows  the  same  general  shape  as 
the  experimental  curve.  Figure  40  shows  the  dependence  of  the  direc¬ 
tional  intensity  elso  on  the  cone  angle  but  for  half  the  axial  length 
used  in  Figure  39.  Fortthis  length,  four  inches,  the  patterns  coincide 
until  the  cutoff  for  the  smallest  cone  angle  is  reached.  Up  to  that 
point,  the  shielding  effect  for  all  three  cones  is  practically  the  same, 
being  at  the  same  time  almost  negligible. 


TABLE  7 

DIFFUSE  LARGE  CONE  PARAMETERS 


Fig. 

No. 

Curve 

No. 

Axial 

Length 

(in.) 

Base 

Rad. 

(in.) 

Open 
Rad. 
(in. ) 

Cone 

Ang. 

(deg.) 

Cone 
Refl . 

A 

Wqi 

(watts/cm.2) 

36 

1 

8 

1.375 

7.919 

39.29° 

.02 

.1  X 

10"3 

- 

2 

8 

1.375 

7.919 

39.29° 

.30 

.1  X 

10-3 

- 

3 

8 

1.375 

7.919 

39.29° 

.60 

.1  X 

10  , 

- 

4 

8 

1.375 

7.919 

39.29° 

.90 

.1  X 

10"3 

37 

1 

2 

1.375 

11.367 

68.35° 

.02 

.1  X 

lO-3 

- 

2 

4 

1.375 

9.992 

68.18° 

.02 

.1  X 

10*3 

38 

1 

8 

1.375 

7.919 

39.29° 

.02 

.1  X 

10~3 

- 

2 

8 

1.375 

7.919 

39.29° 

.02 

.15 

- 

3 

8 

1.375 

7.919 

39.29° 

.02 

.30 

39 

exp.* 

8 

1.375 

2.5 

8° 

.02 

.1  X 

10"3 

- 

1 

8 

1.375 

2.5 

8° 

.02 

.1  X 

10*3 

- 

2 

8 

1.375 

7.919 

39.29° 

.02 

.  1  X 

10"  3 

- 

3 

8 

1.375 

34.8 

76.54° 

.02 

.1  X 

lO*3 

- 

4 

8 

1.375 

81.9 

84.33° 

.02 

.1  X 

lO-3 

40 

1 

4 

1.375 

5.919 

48.66° 

.02 

.1  X 

lO"3 

- 

2 

4 

1.375 

9.992 

68.18° 

.02 

.1  X 

10'3 

- 

3 

4 

1.375 

21.52 

78.77° 

.02 

.1  X 

10"3 

*  Represents  experimentally  determined  curves. 
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Figure  42.  Directional  Intensity  of  Figure  43.  Comparison  of  Aooroxi , 

Diffuse  Double  Cone.  Analyti-  and  Numeri.ai 

Prediction  for  Diffuse 
Cyl  i  ride  r  . 


TABLE  8 


DIFFUSE  CYLINDER  WITH  DIFFUSE  SPHERICAL 
INNERBODY  (FIGURE  41) 


Curve 

No. 

Cyl. 

Length 

(in.) 

Cyl . 
Rad. 
(in.) 

Sph. 

Rad.  R 
(in.) 

Cyl. 
Refl . 

A 

Sph. 

Refl. 

P. 

Pos . 

of 

lb 

W(Cyl.&  Sph.) 
Watts/cm. 2 

1 

7.8125 

2 

1.5 

.02 

.02 

i/2 

.1  x  10“3 

2 

7.8125 

2 

1.5 

.02 

.02 

i-R 

.1  x  1CT3 

exp.* 

7.8125 

2 

1.5 

.02 

.02 

JL-R 

.1  x  10-3 

3** 

7.8125 

2 

- 

.02 

- 

- 

.1  x  10*3 

*  Denotes  experimental  curve 
**  No  innerbody 


Diffuse  Cylinder  with  Diffuse  Spherical  Innerbody 

The  last  family  of  related  curves  is  shown  by  Figure  41.  Curves  1  and  2 
here  depict  the  same  cylinder  and  sphere  but  with  different  positions  of 
the  innerbody.  For  the  particular  values  of  reflectivities  used  here  for 
the  spherical  and  cylindrical  surfaces,  both  being  .02,  the  axial  loca¬ 
tion  of  the  innerbody  has  no  apparent  effect  on  the  shape  of  the  pattern 
since  they  coincide. 

Both  axial  positions  give  the  same  shielding  effect  when  compared  with 
the  same  cylinder  without  an  innerbody.  The  effectiveness  of  the  inner- 
body  is  readily  apparent  in  Figure  41. 

An  isolated  case  is  the  diffuse  double  cone,  whose  experimental  and  pre¬ 
dicted  patterns  are  shown  in  Figure  42. 

EXPERIMENTAL  VS.  PREDICTED  RESULTS 

The  discrepancies  between  the  experimental  and  predicted  patterns  seem  to 
accrue  to  the  experimental  curves.  This  can  be  checked  at  0  =  0  by  cal¬ 
culating  the  normal  radiant  intensity  for  a  given  shield  aperture.  Unac¬ 
counted-for  atmospheric  attenuation  effects  probably  contribute  the  major 
portion  of  the  differences.  Curves  which  support  the  prediction  formula¬ 
tion  for  diffuse  shields  are  shown  in  Figure  43.  The  patterns  predicted 
by  the  numerical  integral  equation  formulation  agree  well  with  the  approxi¬ 
mate  analytic  closed  form  developed  in  Appendix  C. 

On  the  basis  of  the  curves  discussed  above,  the  cylindrical  shield  offers 
the  optimum  shielding  effect  for  a  wider  range  of  aspect  angles  and  for 
higher  values  of  reflectivity. 
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APPENDIX  A 


OBSCURATION  WITHIN  A  CONICAL  SHIELD 
WITH  SPHERICAL  INNERBODY 


Figure  44  shows  the  type  of  shield  configuration  under  consideration, 
namely,  a  truncated  right  circular  cone  with  a  spherical  innerbody  located 
on  the  axis.  For  explicitness,  a  viewing  point  is  taken  on  the  cone  at 
point  x,  and  the  obscuration  from  this  point  is  shown  by  the  shaded  por¬ 
tions  of  the  enclosure.  Attention  will  be  directed  to  the  partially  ob¬ 
scured  ring  elements  on  each  surface.  The  coordinates  used  to  specify 
these  elementary  ring  areas  are  also  shown  in  Figure  44.  (The  dimensions 
of  the  shield  are  shown  in  Figure  46).  Because  of  the  axial  symmetry, 
the  viewing  points  are  chosen  to  lie  in  the  x-z  plane,  allowing  the  extent 
of  obscuration  for  each  ring  to  be  expressed  by  a  single  parameter.  This 
parameter  is  the  z-coordinate  of  the  point  of  intersection  between  the  ring 
element  and  boundary  of  the  obscured  zone.  Such  points  will  be  symbolized 
by  zjj,  denoting  the  obscuration  of  the  jth  surface  as  seen  from  the  ith. 

The  use  of  the  zjj  points  will  be  illustrated  by  an  example.  Consider 
the  partially  obscured  conical  ring  at  f  in  Figure  45,  which  shows  the 
cross  section  through  the  cone  at  ?  parallel  to  the  y-z  plane.  It  is 
seen  from  Figure  45  that  the  azimuthal  angle  is  given  by 

4ll<x,S)^  n.  arccos[ru(x>f)/5(?)i  ,  zu(x,f)<0  (16) 

=  arccos^  Z  ^  (*.f )/«(?}]  .  *n<*.f)>0. 


For  the  conical  ring  at  £  ,  the  ratio  of  the  unobscured  periphery  to  the 
total  is  given  by 


r(x,f)=  )/rr  .  (17) 


The  function  T(x,£  )  very  easily  serves  the  ultimate  purpose  of  this 
analysis,  since  multiplying  the  kernel  for  theunobscured  ring  by  T(x,f  ) 
gives  the  kernel  when  the  ring  is  partially  obscured.  A  similar  function, 
Tr  (x,/0,  serves  th“  same  purpose  for  the  partially  obscured  annular  ring 
on  the  end  disk. 

However,  a  slightly  different  approach  is  used  when  the  sphere  is  viewed 
or  when  the  viewing  point  is  on  the  sphere.  It  is  necessary,  in  these 
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Figure  44.  Conical  Shield  With  Spherical 
Inner  body 
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cases,  to  put  Eq.  (16)  to  the  form 


U 


=  Tf  - 


arccos  (zjj/rj),  Zjj  *0 


=  arccos(zij/rj ),  z^j  ?  0 


(18) 


where  rj  represents  the  appropriate  radius  of  the  elementary  ring  being 
viewed. 


Z 


Figure  45.  Obscuration  of  Conical 
Ring. 


For  example,  when  viewing  the  sphere  from  the  conical  ring  at  x,  the 
values  of  zi3(x,y)  and  r3  =  Rsini"  are  substituted  into  Eq.  (18)  to 
specify  an  azimuthal  angle  The  angle  ^3  is  in  turn  applied 

directly  in  the  kernel  equation  for  radiant  interchange  between  the 
partially  obscured  spherical  ring  and  the  viewing  point.  (See  Appendix  B.) 
A  similar  approach  is  followed  when  the  sphere  is  viewed  from  the  source 
, disk. 

Thus  the  crux  of  the  obscuration  problem  is  to  derive  analytic  expressions 
for  the  points  z^j  as  seen  from  viewing  points  on  each  surface.  Strictly 
within  the  framework  of  analytic  geometry,  the  analysis  is  straightforward 
but  lengthy;  only  the  mair.  features  of  the  approach  and  the  results  will 
be  given  below. 
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Viewing  Point  on  Cone 

(1)  View  of  Cone  -  For  partially  obscured  conical  rings,  the  analysis  cen¬ 
ters  around  the  imaginary  cone  which  has  its  apex  at  the  viewing  point  at 
x  and  is  tangent  to  the  spherical  innerbody  of  radius  R.  Because  this 
tangent  cone  is  thus  circular,  it  is  possible  to  obtain  zjj  in  terms  of 
the  intersection  of  the  familiar  conic  sections  and  the  circle  of  radius 
$  .  According  to  Reference  17,  the  equation  of  the  tangent  cone  from 
the  point  (x,o,r)  (see  Figure  44)  is  given  by 

2 

fx  -  Jx)2  +  r2  -  R2j  |(x  -  Ji)2  +  y2  +  z^=  [(x  -  i1)(?-  Jt])+rz  -  R2]  (19) 

Equation  (19)  can  be  manipulated  to  yield  the  conic  section  which  is  ob¬ 
tained  when  a  plane  parallel  to  the  y-z  plane  is  passed  through  any  point 
in  the  obscured  region  of  the  cone.  The  result  is  given  in  the  standard 
form  for  the  conic  section)  i.e., 


y2 

“52-  + 


(z  ■  h)2  _  : 


where 


_  R|x  -  (  I  ((x  -  ii)2  +  i2  -  R2]^ 
(x  -  lt)2  -  R2 


R2  (» -r  )2 

(x  -  l!)2-  R2 


H  =  r[*?  -•*1  <*  +  f>+-li2  -  R2] 

(x  -  jlj )2  -  R2 - 


If  the  equation  for  the  circle  is  written  as 


y2  +  z2  =  f2  , 
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then  substituting  Eq.  (21)  into  (20)  and  solving  for  z  gives  the  retired 
locus  of  intersection, 

2U(x,  )  =  hb2  -  Hlfb2  (fc2  •  »2  -  Cl  1  .  (22) 

b2  -  a2 

(2)  View  of  End  Disk  -  When  obscuration  of  the  end  disk  occurs,  the  pre¬ 
ceding  discussion  for  the  conical  ring  holds  intact  here,  except  that 
is  set  equal  to  zero  and  (  )  is  replaced  by  ,  the  radius  of  a  partially 
obscured  plane  annulus  of  the  end  disk.  The  points  of  intersection  are 
thus  given  by  Zj^fx,  )*  where  now 

a  R  x  (x  ~  -*l)2  +  r2  -  R2  ^ 

(x  -  JLX)2  -  R2 


b2  = 


(R  x)2 

(x  -ij)2  -  R2 


h  -  ~  *1  *  -  *21 

(x  -  i,)2  -  R2  ' 

(3)  View  of  Sphere  -  The  obscuration  of  the  sphere  as  seen  from  point 
(x,o,r)  on  the  cone  is  depicted  in  Figure  46.  The  plane  marking  the 
boundary  of  the  visible  segment  forms  a  circular  base  for  the  tangent 
cone,  and  it  appears  as  a  chord  (line  C)  in  Figure  46. 


Figure  46.  Obscuration  of  Spherical  Innerbody. 


79 


The  equation  of  line  C  is 


z  =  R  sin(<-  0) 


-  .i5_-.il)  [f-  Ij 


Rcos«-  0)J 


(24) 


where 


tf  = 


arctan 


0  =  arctan(  [  (x  -  ij  )2  +  r2  -  R2  ]  Vr) 


r  =  a  +  xtan0. 


The  parameter  £  denotes  the  x-coordinate  of  any  point  on  the  line  C. 

The  desired  points  of  intersection  ^3  are  found  by  substituting 
f  =  ii  +  Rcosir  into  Eq.  (24 ) |  i.e., 

Zi3(x,y)  =  (R/r)  [  R  +  (ix  -  x )  cos  Y ]  .  (25) 


The  unobscured  spherical  segment  lies  within  angular  limits 


*1  =  0 
*2  =  *  +  0  » 


where  and  0  are  defined  above. 


Viewing  Point  on  End  Disk 

View  of  Cone  -  When  the  viewing  point  is  the  point  (0,0,^)  on  the  end  disk, 
the  roles  of  y°  and  r  are  merely  interchanged  in  Eq.  (23);  and  setting 
x  =  0  makes  that  equation  applicable  for  Z2l(/»x)* 
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View  of  Sphere  -  The  equation  for  line  C  in  this  case  follows  from 
Eq.  (24)  by  putting  x  =  0  and  replacing  r  by  /  .  Consequently,  the 
points  of  Intersection  are  given  by 


z23(/>»r)  =  (R/P  )  [R  +  licosr]  . 


(26) 


The  angular  limits  of  the  unobscured  segment  are  likewise  modified. 


Viewing  Point  on  Sphere 

View  of  Cone  -  When  the  viewing  point  is  on  the  sphere,  the  visible  por¬ 
tion  of  the  enclosure  lies  above  the  plane  tangent  to  that  point.  The 
viewing  point  is  specified  by  angle  x  when  it  lies  in  the  x-z  plane,  as 
shown  in  Figure  46.  There  the  tangent  plane  at  point  p  is  represented  by 
line  L,  whose  equation  is 


*,,(  r,«)  =  -iii— . 

sinr 


(27) 


Equation  (27)  is  clearly  the  locus  of  the  required  intersection  points. 


View  of  End  Disk  -  The  points  of  Intersection  between  partially 
obscured  annuli  of  the  end  disk  follow  Eq.  (27)  by  setting  x  =  0. 

The  foregoing  analysis  can  be  applied  for  the  reversed  cone,  i.e.,  when 
the  wide  end  is  at  the  origin,  simply  by  reversing  the  sign  on  the  cone 
angle  9  and  for  the  cylinder  by  setting  9  equal  to  zero. 


/ 
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APPENDIX  B 


KERNELS  FOR  CONICAL  SHIELD  CONFIGURATIONS 
WITH  SPHERICAL  INNERBOOY 


Diffuse  Shields  -  Except  for  the  double  cone,  all  the  kernels  for  the 
diffuse  shields  considered  in  this  report  follow  from  those  for  the 
truncated  cone,  wide  end  open,  with  spherical  innerbody.  See  Fibres  44 
and  46  for  the  parameters  of  this  configuration.  Therefore,  in  Eqs.  (28) 
through  (34),  the  kernels  for  the  truncated  cone,  small  end  open,  and  the 
cylinder  are  obtained  by  replacing  the  cone  angle  0  by  -  0  and  zero,  res¬ 
pectively.  For  radiant  interchange  involving  the  sphere,  the  keraels  were 
derived  by  integrating  the  fundamental  equations  (Eqs.  (11a)  or  (lib )] with 
Di(«i)  =  D 2(<2).  The  remaining  kernels  were  derived  by  algebraic  techniques 
(angle  factor  algebra).  These  techniques  are  described  and  applied  in  Ref¬ 
erences  3,  5,  6,  7  and  8.  Hence,  little  is  gained  by  including  the  lengthy 
details  of  the  derivations  here. 

Accordingly,  the  kernels  are  listed  for  radiant  energy  transfer  between 
the  following  elementary  ring  elements* 


(1)  Spherical  Ring  to  Conical  Ring 


K3i(r,x)  =  2R2sinf 


2(aib2“«2*>i)-ai2sinr 

!  sin  ^13  \ 

2W(a^2«t)j2  Js  in2/1 

>  aj+biSinX'  cos/8i3  j 

Pa 

4/rR* 


+ 


f-  t.n‘7 

*  \  l  ai+bisinjr 


(28) 


where 


&  =  a  +  xtanO 


4RS(a1a?-bib?sinT)^a13-2a1b12sinV 
4R  (aj2-  b^sinV)*1 


b:  =  2R* 
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tan  @22  1  I .  (30) 

2  Jj 

'  where 

a3  =  ^1  +  Rcos*)(R  +  X^osT) 

=  -  (i^  +  Rcos*')JP  sinT 
a4  =  U\  +  Rcosr)2  +  /*2  +  R2sin2f 
b4  =  -2/)Rsir,y 


2^4^3-a4a3^ 

(a42-b42)3/a 


a4“b4 

a4+b4 
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(4)  Annular  Ring  to  Spherical  Rlnc 
E»y  the  reciprocity  theorem, 


k23  = 


R2sin  t 


(5)  Conical  Ring  to  Conical  Rinc 


Kll(x»?)  =  —  f  i  -  lx~f>secef(x-f)2sec2fr*6  a2+xftan2fr»-a(x+f)tane]} 

2Sl  lu. 


)^sec^&+-4[a^+xf  tan^£.a  (x+£)tan8j  J 


(6)  Conical  Ring  to  Annular  Rlnc 


Kiotx*/*)  — 


x$(x2+  S2-p*)-x2  tane(x2+  S2+P2) 


{  x2(x2+2^24-2)^)+(J2-/32)2  J: 


(7)  Annular  Ring  to  Conical  Rlnc 


By  reciprocity  theorem. 


Si  =  *  K12 


(8)  Conical  Ring  in  First  Section  of  Double  Cone  to  Conical  Ring  In 
Second  Section  -  See  Figure  47  for  parameters  of  double  cone.  The  remain¬ 
ing  kernels  for  the  double  cone  follow  from  the  above  kernels. 


Kll(xl»S)  “ 


(f2-xi)4+2(?2-xi)2(^l2+<f22)+(<5i2-cT22)2 


,  (35) 
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Figure  47.  Geometry  of  Double  Cone. 
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where 


N  =  4i'i^p1f2(^12-*-^2)-t-(^2-x1)  J(^12+^2)(rtan9^atan»')  - 

2  Sl  $2(£ltan*  +^2tan0^  +(?2‘xl)2[3  $1  $2  + 

(?2-x1)(rtan»+-atan^  )]}  -  [($2  -  xx)2  +  +  <$"22  ]  3  . 

<5^  =  r  +  x^t^n 
S  2  =  ®  ■  ?  2tan0 

(9)  Conical  Ring  in  Second  Section  of  Double  Cone  to  Conical  Ring  In 
First  Section  -  By  the  reciprocity  theorem  and  by  replacing  £2  x2 

and  x^  with  %  ^ , 

c' 

Kh(x2»?i)  ~  — (?l»x2^* 


where 


S1  =  r  +  ^tanV* 

I 

<5  2  =  a  -  x2tan0. 


Specular  Cylinder 

Annular  Diffuse  Ring  of  End  Disk  to  Specular  Cylindrical  Rinc 


[r’,x)=  _4_K(r’,xn: 

L  An 


where 


_  (2n-l)a+r’ 

- ; -  -  A 

a+r  * 
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APPENDIX  C 


APPROXIMATE  ANALYTIC  FORKJLATION 
FOR DIRECTIONAL  RADIANT  INTENSITY  OF  A  DIFFUSE  CYLINDER 


The  approximate  analytic  solution  given  by  Moon  for  the  diffuse  circular 
cylinder  is  of  the  form 

Wj(x)  =  B  cosh  kx  +  Csinh  kx  +  D 


W2(r* )  =  const. , 


(38) 


where  W^(x)  is  the  apparent  emittance  distribution  as  a  function  of  length 
x  along  the  cylinder  wallj  B,  C,  and  D,  are  constant  and  depend  on  the  shield 
parameters  -  length,  radius  reflectivity,  and  initial  emittance  and 
*2(0  is  the  final  emittance,  including  interreflections,  which  is  assumed 
constant  over  the  source  disk.  (The  mathematical  equations  for  B,  C,  D, 
and  W2  are  somewhat  lengthy,  and  their  reproduction  here  would  add  nothing 
essential  1  see  Reference  12,  Eqs.  (21  through  (24).) 

The  coordinate  system  of  the  cylindrical  shield  is  shown  in  Figure  29.  It 
is  easily  seen  from  Figure  29  that  source  disk  is  not  visible  for  aspect 
angles  greater  than  0C  =  arctan(2a/fc ).  Figure  48  shows  the  interior  of 
the  cylindrical  enclosure  projected  onto  a  plane  perpendicular  to  the 
line  of  sight  at  angle  0.  The  rectangular  coordinate  system  s,t  lies  in 
this  plane.  . 


Figure  48.  Projection  Plane  Coordinates. 
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Consider  first  the  c«se  0>0C.  The  essence  of  the  problem  is  to  relate 
the  emittanee  function  Wj^(x)  in  the  shield  coordinate  system  to  a  func¬ 
tion  L(s,t)  in  projection  plane  coordinates  and  then  to  integrate  the 
function  L(f,t)  over  projected  area  Ap  of  the  cylinder.  In  short,  the 
radiant  intensity  at  angle  0  will  take  the  form 


L(s,t)dA  . 


(39) 


The  function  L(s,t)  is  determined  by  the  following  arguments  on  the  cylin¬ 
der  the  emittanee  is  constant  around  each  elementary  ring,  which  shows  up 
as  an  arc  of  an  ellipse  on  the  projected  area.  Therefore,  the  point  (s,t) 
can  be  considered  to  lie  on  an  elliptical  arc,  the  equation  of  the  ellipse 
being 


i»-h)2  +_tL=  1. 

(acos0)2  a2 


(40) 


From  Eq.  (40)  the  distance  h  is  given  by 


h(s,t) 


-  acos0  |^1  -  _JL_ 


(41) 


Since  the  emittanee  is  constant  along  the  parent  circular  ring,  it  is 
also  constant  along  the  corresponding  elliptical  arc.  Hence,  it  is 
possible  to  relate  the  function  L(s,t)  to  its  value  at  the  point  (so,0), 
where  sQ  is  the  point  of  intersection  of  the  elliptical  arc  and  the 
s-axls.  By  symmetry,  this  point  is  in  the  horizontal  diametral  plane  of 
the  cylinder.  From  Figure  48,  sQ  is  equal  to  h  -  acos0  and,  with  Eq.  (40), 
gives 


s0  =  s  -  acos0  [  1  -  (1  -  t2/a2)^  ] . 


(42) 


By  defining 


g(t)  =  acos0  [l  -  (1  -  t 2/*2)t  J 


(43) 
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it  is  possible  to  set  up  the  following  functional  relationship! 


L(s,t)^  L(sq,0)  X  Lfs  -  g(t )} 


(44) 


Before  Eq.  (44)  can  be  useful,  a  functional  relationship  must  be  estab¬ 
lished  between  th*s  coordinates  of  the  (s,t)  system  and  the  x-axis  of  the 
cylinder.  From  Figure  48*  it  is  easy  to  establish  the  following  corre¬ 
spondence. 


L(0,0)  ^WiU  -  —2 - ) 

tantf 


L(s,0)^W,(-  — § —  +  £  - 

sin0 


LCs.t)^*!  f-  .ilgiLl,  +  X 

l  sin0 


(45a) 

(45b) 

(45c) 


Through  Eqs.  (45c)  and  (43),  it  is  possible  to  transform  the  emittance 
function  along  the  x-axis  of  the  cylinder  to  a  distribution  function  along 
the  x-axis  of  the  projected  area,  and  ultimately  to  the  entire  projected 
area. 

But  it  is  necessary  yet  to  transform  the  cylindrical  surface  element  dA 
into  a  corresponding  differential  projected  area  element  dff*  .  In  terms 
of  the  cylindrical  coordinates  of  Figure  29,  dA  is  given  by  ad9dx.  From 
Figure  48,  it  is  found  that  dx  =  ds/sin0  and  ad9  =  dt/cos0j  hence, 

dA  =  adxdG  =  —S.—t _ =  ~6<r  .  (46) 

sin0cos0  sin0cos© 


At  this  point  it  is  pertinent  to  consider  the  angle  <x  made  by  the  line  of 
sight  and  the  normal  to  any  element  of  area  on  the  inner  cylinder  wall. 

For  the  diffuse  cylinder,  the  radiant  intensity  from  each  differential 
area  element  in  the  direction  of  observation  is  given  by 

dJ(«)  =  dJn(*)C0.«  =  -iliiiil  -C°f«  d'  (47) 

^  sin0cos© 
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where  dJn(x)  represents  the  intensity  along  the  normal  to  each  differen¬ 
tial  element.  It  was  possible  in  the  laboratory  experiments  to  consider 
the  lines  of  sight  to  lie  in  planes  parallel  to  the  diametral  plane  of 
the  cylinder.  Under  this  condition,  the  direction  cosines  of  the  line 
of  sight  are 


Ul  =  sin0,  Jt2  =  0,  J3  =  -  cos0). 


whereas  those  of  the  normals  to  the  inner  cylinder  wall  are 


(mi  =  cosd,  m2  =  -  sin©,  17)3  =  0) | 


hence,  the  cosine  of  the  angle  between  the  line  of  sight  and  the  normals 
is 


cos*  =  mi  li  +  ^2^-2  +  m3^3  =  cos©  . 


(48) 


Equation  (48)  leads  to  an  interesting  result;  for  when  cos*  as  expressed 
by  Eq.  (48)  is  substituted  into  Eq.  (47),  it  is  seen  that  the  angular 
dependence  of  the  radiant  intensity  from  the  entire  projected  area  is  due 
only  to  the  aspect  angle  0.  With  the  results  of  Eqs.  (45c),  (46),  (47) 
and  (48),  Eq.  (39)  is  transformed  into  the  following  prediction  formula: 


s~9(t)  +  i 

s  in0 


-  r  Ids  dt, 
tan0J 


(49) 


where  it  is  understood  that  the  variable  x  of  Eq.  (38)  is  to  be  replaced 
by  the  argument  within  the  braces  of  Eq.  (49)  and  that  the  integration  is 
then  to  be  taken  over  the  projected  area. 

The  problem  generated  by  substituting  Eq.  (38)  into  the  prediction  formula 
is  somewhat  complicated;  but  when  the  cylinder  is  maintained  at  a  sub¬ 
stantially  lower  temperature  than  the  source,  Eq.  (38)  reduces  to  a  simple 
exponential,  i.  e., 

"Ki.b)(x)  =  Be'kx  •  <50) 


91 


In  the  form  required  for  Eq.  (49),  Eq.  (50)  is  given  by 

W1(iab)(s*t)  =  B  expf  -k[j-(a2-t2)^  sec0  -  s  csc0]J. 


The  double  integration  of  Wj.  (iafc)  (s  »t)  is  tedius  but  straightforward,  and 
the  details  are  omitted.  The  final  result  is  given  by  the  following: 


J(0)  = 


_  2  Be 


~[f 


exp  »Ca-*2)*  dt-a  1 , 
tanG  I 


0>0C  . 


There  is  no  essential  difference  in  the  derivation  when  part  of  the  source 
disk  is  visible;  i.e.,  the  final  prediction  formula  for  this  case  is 

J(0)  =  2Be  klsinG-  [a(pekl-n  ♦ 
k 


2k(a2-t2)* 

tan0 


dt  +  *2^(0), 


where 


■(" 


tan20 

tan^ 


and  0  <  0 
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